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[1] Teaching Mathematics with DERIVE and the 77-92 , Proceedings of the International 
DERIVE and TI-92 Conference 1996 at Schloss Birlinghoven, Bonn. 

543 pages with 54 contributions (13 German, 39 English), available at ZKL, Prinzipalmarkt 38, 
D 48143 Muenster, Germany; email: bemtz@uni-muenster.de. (DEM 50.00) 

[2] La calcolatrici grafiche nelPinsegnamento dell matematica, G.C.Barozzi, S Cappuccio 
146 pages, Pitagora Editrice Bologna 1997, ISBN 88-371-0891-5 

This new booklet is addressed to Secondary School teachers and deals with graphic and 
programmable calculators (most TI-92) in Maths teaching and learning. 

(Arithmetic, Algebraicallanguages, Calculators as a tool for Algebra, C. as an aid for learning, 
The discovery method, C. in problem solving, C. as a blackboard, Geometry, Locus of points: 
conics, Programming, Lists, Geometric transformations. A "Matura" problem solved via TI-92) 

[3] Mathematical Activities with DERIVE, Edited by E. Graham, J.S. Berry and A.J.P. Watkins, 
207 pages, Chartwell-Bratt 1997, ISBN 0-86238-478-8 

This book brings together a range of interesting examples of the use of DERIVE which have 
been contributed by a number of DERIVE users (22) from around the world. 

(What was the question?, DERIVE and Mechanics, Linear Algebra Using DERIVE and 
Methods of Assessment, A Modelling Workshop, The Zeros of a Cubic Polynomial, Examples 
of Studies of Simple Dynamic Systems, Construction of Fractal Sets, Concepts of Differential 
and Integral Calculus and others. ) 

[4] The State of Computer Algebra in Mathematics Education, Edited by I.S. Berry, 

J. Monaghan, M. Kronfellner, B. Kutzler, 213 pages, Chartwell-Bratt 1997, 

ISBN 0-86238-430-3 

A report of the International Symposium on Computer Algebra in Mathematics Education held 
in Hawaii, August 1995. 

(The Curriculum, Assessment, The design of CAS written support material, Attracting Teachers 
to Use CAS in the Classroom, Understanding) 

[5] A Survey of Models for Tumor-Immune System Dynamics, Editors J.A. Adam, N. Bellomo, 
344 pages, Birkhauser Boston 1997, ISBN 0-8176-3901-2 

This unique book is a collection of seven interdisciplinary surveys on modeling tumor dynamics 
and interactions between tumors and the immune system. It is a excellent resource and survey 

for applied mathematicians, mathematical biologist and biologists interested in modeling 
methods in immunology and related sciences. 

Among others you can find various Diffusion Models, The Predator-Prey Approach, A Model 
of Tumor Cell/Immune System Interaction, Comments on Catastrophe Theory , Mathematical 
Modeling of Tumor Gromh Kinetics, Tumor Gromh as a Dynamical System, Logistic Model, 
Multicell Spheroid Model. References for further reading 


Interesting WEB sites 


http://www.imag. fr/imag/cabri.html 
http://www.ti. com/calc/docs/cabri .htm 

http://forum.swarthmore.edu/ 
http://www .cut-the-knot.com/ 


from the Cabri developers in Grenoble 
from Texas Instruments 

including The Math Forum, Ask Dr. Math, 
Internet Resources, Cabri Information, ... 

(cti, Birmingham) 

Interactive Mathematics and Puzzles by 
Alexander Bogomolny. An award winning Math 
site on a variety of topics, recommended by Jan 
Vermeylen 


Please let us know your favourite site! 
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Liebe DERIVE- und Tl-Freunde, 


Dear DERIVE- and Tl-friends, 


Der letzte Schritt in der Fertigstellung eines je- 
den DNL ist immer das Verfassen des Letters. Bis 
zum letzten Augenblick erreichen mich interessante 
Nachrichten und Noor, meine liebe Frau, drangt und 
drangt: Jetzt mach doch endlich fertig!" 

Aus den letzten drei Monaten habe ich beson- 
ders schone Erinnerungen an die gelungene DERIVE- 
und Tl-Tagung in Munster. Herzlichen Dank von die- 
ser Stelle an die Organisatoren mit Dr. Berntzen und 
Barbel Barzel an der Spitze. Auch T A 3-Europe nimmt 
schon langsam Form und Gestalt an. Wenn Sie noch 
nicht wissen, was T A 3 ist, dann besuchen Sie bitte die 
Homepage von Bert Waits. 

Ich freue mich aber auch schon 
auf das vor uns liegende Treffen in 
Schweden Die FUN-Conference (Fun 
in Teaching Mathematics) versammelt 
vom 7. bis 9. August mehr als 30 Per- 
sonen aus vier Kontinenten zu einem 
Austausch ihrer Ideen und Arbeiten. 

Viele Freunde werden sich wieder 
treffen und dort sicherlich neue 
Freunde gewinnen. Ich finde es sehr 
wichtig, dass, sich der Kreis immer 
wieder erweitert und erneuert. Damit 
ist doch sicher gestellt, dass wir nicht 
in unserer Routine erstarren. Da auch 
David Stoutemyer und Bert Waits 
dabei sind, konnen wir sicher sein, 
dass unsere Vorstellungen sofort an 
die richtige Adresse gelangen. 

Vor der FUN-Konferenz gibt es am selben Ort 
(Kungsbacka) eine von David Sjostrand gemeinsam 
mit der DUG organisierte Lehrerfortbildungsveranstal- 
tung, bei der 50-60 Teilnehmer erwartet werden. 

Und auch die „gro3e“ Internationale DERIVE und 
TI-92 Konferenz 1998 in Gettysburg wirft schon ihre 
Schatten voraus. Beachten Sie bitte die Information 
ab Seite 1 1. Car! Lein bach hat sich bereits viel Mu he 
gemacht, ein wunderbares Nebenprogramm zusam- 
men zu stellen. Mit viel Freude habe ich von ihm ge- 
hort, dass, er nach seinem Herzleiden wieder voll 
hergestellt ist. Carl und Pat, wir wunschen Euch alles 
Gute! Fur den fach lichen Teil sind wir alle verantwort- 
lich. Ich lade Sie herzlich ein, an dieser Konferenz als 
Vortragender, als Leiter eines Workshops oder ein- 
fach nur als Zuhorer teilzunehmen. 


The last step completing each DNL is writing 
the Editor's Letter. Interesting news reach me until 
the last moment and my wife Noor is urging and 
urging. "It's time to come to an end!" 

I have best memories on the excellent DERIVE 
and TI-92 Conference in Munster, Germany. Many 
thanks from here to the organizers led by Dr. 
Berntzen and Barbel Barzel. T A 3 -Europe is now 
beginning to show shape and form. If you don't know 
what T A 3 is then please visit Bert Waits' homepage. 

I am looking forward to the meeting in Sweden 
ahead of us. The FUN-Conference (Fun in Teaching 
Mathematics) will gather more 
than 30 persons from four 
continents to exchange their ideas 
and work. Many friends will meet 
again and will make new ones. 
From my point of view it is very 
important to enlarge and to renew 
the circle. That will ensure that 
we will not stiffen in our routine. 
David Stoutemyer and Bert Waits 
will be with us, so we can be sure 
that our ideas and proposals will 
immediately reach the right 
address. 

Just before the FUN- 
Conference a Teachers' Training 
Conference will take place in the 
same location (Kungsbacka) organized by David 
Sjostrand and the DUG. We are expecting 50 -60 
participants. 

And above all, the "big" International DERIVE 
and Tl-Conference casts its shadow ahead of it. 
Please note the information on page 11. Carl 
Leinbach has spared no effort to arrange a wonderful 
Spouses’ Program. I was very glad to hear from Carl 
that he has fully recovered from a severe heart attack. 
Carl and Pat, the best to you from all of us. We all are 
responsible for the technical part of the Conference. 
So I'd like to invite you all to attend to give lectures 
or to lead workshops or to participate and help to 
make this Conference as successful as the former 
ones have been. 



Allen Freunden mit Internet-Anschluss, mochte 
ich die angebotenen Webseiten warmstens empfeh- 
ien. Sie entpuppen sich als wahre Fundgruben. 

Ich wunsche Ihnen alien einen schonen Sommer 
mit erholsamen Ferien - unseren Freunden auf der 
anderen Halbkugel einen nicht zu unangenehmen 
Winter - und freue mich bereits auf den nachsten DNL 
mit hoffentlich vielen interessanten Neuigkeiten aus 
der Welt von DERIVE und TI-92. 


Josef 



I warmly recommend the websites presented on 
the information page. You will find great resources 
and many links to other sites. 

I wish you all a wonderful summer with 
relaxing holidays - our friends living in the other 
hemisphere I wish a not too bad winter - and I am 
looking forward to the next DNL which will bring 
again many interesting news from the world of 
DERIVE and the TI-92. 

Josef 





EDITORIAL 
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The DERIVE-NEWSLETTER is the Contributions: 

Bulletin of the DERIVE User Group. It is Please send all contributions to the Editor, 

published at least four times a year with Non-English speakers are encouraged to 

contents of 40 pages minimum. The goals write their contributions in English to 

of the DNL are to enable the exchange of reinforce the international touch of the 

experiences made with DERIVE as well as DNL. It must be said, though, that non- 

to create a group to discuss the possibilities English articles will be warmly welcomed 

of new methodical and didactical manners nonetheless. Your contributions will be 

in teaching mathematics. edited but not assessed. By submitting 

We include now a section dealing with the articles the author gives his consent for 

use of the TI-92 and we try to combine reprinting it in DNL. The more 

these modem technologies. contributions you will send the more lively 

■ i and richer in contents the DERIVE 

Editor: Mag. Josef Bohm Newsletter will be. 

A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone: 43-(0)660 31 36 365 
e-mail: nojo.boehm@pgv.at 


Preview: Contributions for the next issues 

3D-Geometry, Reichel, AUT 
Algebra at A-Level, Goldstein, UK 

Graphic Integration, Linear Programming, Various Projections. Bohm, AUT 
A Utility file for complex dynamic systems, Lechner, AUT 
Examples for Statistics, Roeloffs, NL 

Solving Word problems (Textaufgaben) with DERIVE, Bohm, AUT 

Line Searching with DERIVE, Collie, UK 

About the "Cesaro Glove-Osculant", Halprin, AUS 

Hidden lines, Weller, GER 

Fractals and other Graphics, Koth, AUT 

Experimenting with GRAM-SCHMIDT, Schonefeld, USA 

Implicit Multivalue Bivariate Function 3D Plots, Biryukov, RUS 

The TI-92 Section, Waits a.o. 
and 

Setif, FRA; Vermeylen, BEL; Leinbach, USA; Halprin, AUS; Speck, NZL; 
Weth, GER; Wiesenbauer, AUT; Aue, GER; Propper, GER; Koller, AUT; 
Stahl, USA; Mitic, UK; Tortosa, ESP; Santonja, ESP; Wadsack, AUT; 
Schorn, GER, and 

Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag.Josef Bohm 
Herstellung: Selbstverlag 
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Scott A. Guth, Mt San Antonio College 

Hello all! 

I know that this has come up recently, but I’ve lost the response to the question which was asked. 
To keep it simple, how can I fix the expression 

[ random (10), random (10)] 


so that the same numer is returned in both positions of the vector? 


pstoffel @ uia.ua.ac.be 

Albert D. Rich suggested the following method about a year ago (Tue 14 May 1996 for the 
mathematical correctness): 

f (x) := [ x , x ] 

then simplify 

f (random (10) ) 

and then DERIVE returns the desired answer! 


#1: 

[RANDOM (10), RANDOM (10)] = 

[3, 

5] 

#2: 

[RANDOM (10), RANDOM (10)] = 

[2 , 

5] 


fCO := [x, x] 



#4: 

f (RANDOM (10)) = [5, 5] 



#5: 

f (RANDOM (10)) = [8, 8] 



#6: 

■f (RANDOM (10)) = [4, 4lHM 




IFI 7 T 1 Y F£T Y FjT Y FHt Y FE Y Ffi t T "| 

| t |fl 1 gebra |c a 1 c lother |Pr gm 1 0 |c 1 ean Up | 

■ Oand( 18) rand(10)> 

CS 2> 

■ tx x> ■* f(x) 

□one 

■ f(randC10)) 

€5 5> 

■ f(randClG)) 

€6 61 

f <randaG>>l 

liDYREG RAD liUTD 

FUNC H/30 BHTT 


[randInt(l,lo) randlnt/l, 1C))] 

[7 10] n 

/M : =[* x] 

Per tig 

/(randlnt(ljio)) 

[4 4] 

/(randlnt(lj 1 o)) 

[1 1 ] 

//randlntfl, 10j) 

[lO 10] 

1 


DNL: As you can see the same procedure is necessary working on the Tl-handhelds and 
with TI-NspireCAS as well. 
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.Tan Vermevlen, Kapellen, Belgium 

1) While talking in class about singular conic sections I was embarrassed to see that y 2 = 4 is plotted 
with DfW 4.02 as two VERTICAL lines. However, y = 2 is correctly plotted in horizontal 
direction. This error does not occur in the DOS-version. 

2) The expression TERMS (x A 2 + 5y A 2 + i x) is simplified to one term, and I hoped to get 
three of them. 

3) More than once I get remarks from teachers that the soLve -function of equations should have 
the possibility to include/exclude complex solutions. With younger pupils, this is necessary. You 
can make this choice in the Factor-function, why not in soLve? More in general, why is a 
domain declaration of a variable not taken into account while solving an equation in that variable 
or plotting an expression in that variable? 


PA/L: 

1 ) DfW 4.05 plots y 2 = 4 as two vertical lines but DfW 4.04 

lines. 

(In the meanwhile DERIVE 6 plots two horizontal lines.) 

2) In order to decomposite the expression x 2 + / + i-x into 
the variables as Complex. This is the trick 
(This is not necessary in DERIVE 6.) 

2 2 r 2 2 1 

#2: TERMS (x + 5-y + l-x) = Lx , 5-y , i.-xj 

Jan Vermevlen, Kapellen, Belgium 

I give you here some extracts of my correspondence with Theresa Shelby from SWHH about plotting 
problems with DERIVE for Windows. 

Jan: Dear Theresa, 

Allow me to mention a problem I have when plotting a vector of expressions: 

Unlike it is said in the Online Help one of the expressions is plotted in the same color as the 
background. So it is invisible! You can check this with VECTOR (y=k, k, -4, 4, 1/4). 

Theresa: On my computer this does not happen. Some of the colors are very, very close to 
the background color, however. We have a deficiency in our algorithm that only rejects a 
color if exactly the same as the background; we need to add some kind of buffer. There 
could also be a difference in what color value is and what the device actually displays. Any 
way, we do need to work on this area of 2D plotting. Thanks for letting us know. 

Jan: Can I add to the plot problem y 2 = 9 which I reported last week a similar one? 

(x + y - l) 2 = 0 is not plotted at all! 

Theresa: Yes, I tried this also on DfD and only a small segment is plotted and only when the 
scale is 2 : 2. As you noticed it is not plotted at all in DfW. 

I have more information on the last plotting problem you reported - this answer comes via 
David Stoutemyer ... 


recognizes y 2 = 4 as two vertical 


three terms you have to declare 
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This is a double root line. It is a parabolic valley z = (x + y - if that barely touches z = 0, and 
none of our sample points happen to occur where z = 0. It is an unavoidable limitation of 
plotting contours from a finite number of samples. Try (x + y - if or (x + y - if. 


Jan: My remarks: 

The y 2 = 9 problem is not yet solved or answered. I can understand the (x + y - l) 2 = 0 problem but 
when teacing conic sections both problems give a ridiculous impression to the pupils. 

Some days later . . . 

Jan: While looking back at them another remarkable thing came up. 
y = 3 plots one horizontal line as it is supposed to do 


BUT 

y 1 = 9 or y A2 -9 = 0 plot wrongly two vertical lines 


BUT 

0 * x + y 2 = 0 plots correctly two horizontal lines! 

DNL: So does DERIVE for DOS 4. 


Troels Ring, Aalborg, Denmark 

Can anyone explain how to use DERIVE for exploring the general relationship between the correlation 
between two variables and the correlation between functions of two variables. The real problem is to 
see if the correlation between log(A') and log( Y) is a simple function of the correlation between X and 
Y. 

Thanks in advance. 


Julio Gonzales Cabillon, Urugav 

We would like to inform you that the book titled 

“Italian Research in Mathematics Education” 

is now available on-line at our mathematics education server. It compiles the Italian papers presented 
at ICME8, Sevilla. 

URL: http :// 157.253.25.2/ servidor/ em/recinf/libros .html 


ctimath @ bham.ac.uk 

The Third Edition (March 1997) of the MathSkills Newsletter is now freely available in the Newsletter 
archive at the MathSills website: 

URL: http://www.hull.ac.uk/mathskills/newsletters(index.html 


(The link is still valid!) 
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warrensl @ accucomm.net 

Is there a way to create log-log graphs using DERIVEl 

DNL: Try this: 


t is a 2 columns-table ■+ loglog-table oft 

loglog(t) := VECrOR(|"LN(t ), LN(t )], i p DIM(t)) 

#1: L i,l i,2 J 

f is a given function ■+ loglog-transformation of this function 

loglogf(f) := SOLVE( Him Him (LN(LHS(f)) = LN(RHS(f ))) , y) 
#2 : y-?e A v xh^ a j 


1.57 

#3: tab := TABLE (4.2 -x , x 5 50, 400, 50) 


50 

1952.685031 

100 

5797.613911 

150 

1.095752744-10 

200 

4 

1.7213389-10 

250 

2.443515532-10 

300 

3.253341553-10 


#5: log log (tab) 


3.912023005 

7. 

576960643 

4.605170185 

8. 

665201717 

5.010635294 

9. 

301781936 

5.298317366 

9. 

75344279 

5.521460917 

10 

.10377816 

5.703782474 

10 

.39002301 

5.857933154 

10 

.63203957 

5.991464547 

10 

.84168386 


4 

350 4.144156235-10 

4 

. 400 5.110736343-10 


1.57 

#7: loglogf(y = 4.2-x ) 


= 1.57- li + 1. 435 0845 25 




4 


6 


S 


10 
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Scott Guth, USA 

After applying some simple algebra to some trite phrases and clinches as a new understanding can be 
reached of the secret to wealth and success. 

Here it goes. 

Knowledge is Power 
Time is Money and, 
as every engineer knows, 

Power is Work over Time. 


So, substituting algebraic equations for these time-worn bits of wisdom, we get: 


( 1 ) 

K = P 

( 2 ) 

T = M 

( 3 ) 

W 

P = — 

T 


Now, do a few simple substitutions: 

Put WIT in for P in equation (1), which yields: (4) 

Put M in for T into equation (4), which yields: (5) 

Now we’ve got something. Expanding back into English, we get: 

Knowledge equals Work over Money 

What this MEANS is that: 


vv 

K = — 

T 


K = 


W_ 

M 


1. The More you Know, the More Work You Do, and 

2. The More you Know, the Less Money You Make. 

W 

Solving for Money, we get: (6) M = — 


Money equals Work over Knowledge 

From equation (6) we see that Money approaches infinity as Knowledge approaches 0, regardless of 
the Work done. What this MEANS is: 

The More you Make, the Less you know. 

Solving for Work, we get: (7) W = M • K 

Work equals Money time Knowledge 

From equation (7) we see that Work approaches 0 as Knowledge approaches 0. What this MEANS is: 

The stupid rich do little or no work 


Working out the socioeconomic implications of this breakthrough is left as an exercise for the reader, 
Tosh 


p 
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Should Computer Algebra Programs Use 
In xor In |x| as their Antiderivative of 1/x? 

David R. Stoutemyer, SWHH, Hawaii, USA 


The DERIVE™ computer algebra system currently uses In xas its antiderivative of 1/x. 
About once a month we receive a suggestion that we should use In |x| instead. Many current 
calculus texts alternate between these alternatives with no explanation of why. This article is 
intended to provoke sufficient debate so that implementors can decide what to do about this. 


1. Logarithms in Computer Algebra Systems 

For reasons of practicality, most hardware and software regards an expression such as 
ln(z) as denoting only one branch of the infinitely-branched logarithm function. With a 
computer algebra system, users can always specify an infinitely-branched logarithm by 
entering expressions such as ln(z) + 2 n-u /, where n denots an indeterminate integer. After 
dealing with the consequent clutter they will probably decide that this is not the panacea they 
hoped. 

If In denotes one branch, correct use is more likely if In consistently denotes a 
particular branch, such as the principal branch usually defined by 

Re(ln z) = In |z| and lm(ln z) = phase(z), where -n < phase(z) < n. 

Among other prohibitions, consistency to this principal branch entails resisiting the 
temptation to “simplify” ln(exp w) to w unless the algebra system can deduce that 
-n < lm(w) < n. 

2. Why In z is the better choice for complex z 

Most computer algebra programs were not written solely to meet the needs of 
beginning real-variable calculus. It is a desirable goal to have integration give correct results 
over contours in the complex plane as well as over real intervals throughout which the 
integrand is well behaved. The multiple-branched logarithm is analytic except at z = 0, and 
the principal branch is applicable to contours that also don’t cross the negative x-axis without 
returning. 

Consequently, integration of 1/x from a point such as -3 to point such as +2 would yield 
In 2 - (In 3 + /■ n) = ln(2/3) - /'■ n. This is consistent with taking a contour from x = -3 to x = -r 
with r> 0, proceeding along a semicircle of radius r from angle n to angle 0, then proceeding 
along the positive x-axis to x = 2. There is no need to take the limit as r approaches 0. 

, x , d , , | sign(z) 1 , . . , , , 

In contrast , — In z = = where z= sign(z) • |z|. 

dz z conj(z) 


z = conj(z) only for real z, so ln|z| is generally wrong for complex z. 
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3. Why In x is correct even for negative real intervals 

Integrating a continuous real function over a real interval should give a real result. 
However, this does not require that an antiderivative used to determine such an integral must 
be real: Valid antiderivatives can differ from each other by an arbitrary constant, and that 
constant can be complex, such as i-n. Any such constants cancel out in a definite integral, so 
using In x gives the desired real result when integrating between two negative values. 


4. Why In |x| is wrong even for real x 

Here is an excerpt from "A First Course in Calculus" by Serge Lang (Addison-Wesley, 
pp. 136-137): 


"We agree throughout that indefinite integrals are defined only over intervals. Thus 
in considering the function 1/x, we have to consider separately the cases x > 0 and x 
< 0. For x > 0, we have already remarked that log x is an indefinite integral. It turns 
out that for the interval x < 0 we can also find an indefinite integral, and in fact we 
have for x > 0, 


dx = log(-x) . 


Observe that when x < 0, -x is positive, and thus log(-x) is meaningful. The fact that 
the derivative of log(-x) is equal to 1/x is true by the chain rule. 

For x < 0, any other indefinite integral is given by 

log (-x) + C, 

where C is a constant. 

It is sometimes stated that in all cases, 


j — dx = log |x| + C . 


With our conventions, we do not attribute any meaning to this, because our 
functions are not defined over intervals (the missing point 0 prevents this). In any 
case, the formula would be false. Indeed, for x < 0 we have 


j — dx = log |x| + C x 


and for x > 0 we have 


[ — dx = log |x| + C 2 
J x 


in all cases. 





p 10 


David R. Stoutemyer: In xor ln|x|? 


D 


— N— L#26 


We prefer to stick to our convention that integrals are defined only over intervals. 

When we deal with the log, it is to be understood that we deal only with the case 
x > 0." 

Note also that the integral of 1/x diverges for contours that go through x = 0. If users 
are thinking of such a path, then the complex result of using In x will probably alert them to 
their misuse whereas the incorrect real result of using In |x| is more likely to go unnoticed. 

It is true that using In |x| fortuitously gives the correct Cauchy principal value for all real 
intervals. However, Cauchy principal values are customarily and wisely not taught in elemen- 
tary calculus. Their valid use generally requires supporting physical justification, such as the 
fact that the singularities are known to be artificial canceling artifacts of the particular 
mathematical model of a physical problem. 


5. Why some people might nonetheless prefer In |xj: 

Probably the reason some prefer In |x| is aesthetic: For reasons of parsimony we wish 
that real-variable calculus could be done entirely without excursions into complex variables. 

I doubt that this is possible for sufficiently complicated examples. Moreover, even for the 
simple example of integrating 1/x, using In |x| to avoid complex variables unnecessarily 
teaches bad habits. 

Calculus students are introduced to complex numbers in high-school algebra, and they 
are taught that antiderivatives can differ by arbitrary constants. Why not exploit logarithms of 
negative numbers to reinforce both ideas? This is not to say that we should generalize all 
elementary calculus to complex variables -- only that if a complex number is most correct 
and natural in a few places, exploit it rather than fight it. 

It is true that the i-n from In x does not disappear when instead of evaluating a definite 
integral on a negative interval we are merely substituting a single negative number into this 
antiderivative. However, since antiderivatives can differ by a constant, it is hard to imagine a 
good reason for only substituting one number into an antiderivative -- one might as well pick 
a random number. 

On the other hand, there are good reasons to plot an antiderivative as its variable 
varies over some interval, and most 2-dimensional plot packages intended for real values 
conservatively plot nothing wherever a value is unreal. However, it is always possible to plot 
both the real and imaginary parts of an antiderivative such as In x. Indeed, doing so has the 
educational value of helping to reveal the phase discontinuity across x= 0. 


6. Why In |xj is abhorred by computer-algebra implementors even for x positive or 
for x negative: 

It is hard to automate powerful correct automatic simplification of absolute values. 
Once an expression is contaminated by them, they tend to persist unnecessarily in most 
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subsequent derived expressions. For example, many computer-algebra systems don’t 
simplify the real derivative of In |x| all the way to 1/xfor real x. Absolute values also tend to 
thwart operations such as equation solving or further integration. For example, not many 
systems can automatically solve 2|3x — b| > 5 or compute antiderivatives for 1/|x| or 
ATAN |ax+ b|. 


7. Conclusions 

Despite the above arguments, for those users who insist on including an absolute 
value we are considering having a future version of DERIVE include a switch: When In |x| is 
selected, DERIVE would use In \u\ in antiderivatives whenever DERIVE can deduce that u is 
real. Otherwise DERIVE would use In x. Should we include such a switch? If so, what should 
be the default setting? 


As you can see on the following screen shots, the TI-92 offers the "switch" David spoke 
about in the last paragraph. 




Josef 
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Request for Submissions 

Third International DERIVE and TI-92 Conference 
July 14-17 on the Campus of Gettysburg College 
Gettysburg, Pennsylvania 


Submissions are requested for the Third International DERIVE and 77-92 Conference to be 
hedl July 14-17 on the campus of Gettysburg College in Gettysburg, Pennsylvania, USA. 
Submissions may be in form of proposals for 25 or 50 minute presented papers, or proposals 
for 2 hour workshops. All proposals will be refereed by the Conference Co mm ittee who will 
select from among the submissions. Acceptance of submissions will be based on the reports 
of the referees and the availability of program slots. In past conferences there has been an 
approximately 70% acceptance rate of proposals. 

Topics for presented papers are expected to include discussions of at least one of the follow- 
ing topics: the use of DERIVE and/or the TI-92 in the classroom with an analysis of its use; 
applications of DERIVE and/or the TI-92 to topics in school and university level mathematics; 
discussions of the educational issues involved in using a Computer Algebra System (CAS) in 
the teaching of mathematics; the effects of using a CAS on the mathematics curriculum; or 
suggestions for the future design of a CAS for use in the teaching of mathematics. Other top- 
ics for papers may be considered if they are deemed appropriately related to the themes listed 
above by the Conference Committee. 

Workshops may be either introductory tutorial, advanced tutorial, or on the use of DERIVE 
and/or the TI-92 in the teaching of a subject or group of subjects found in the mathematics 
curriculum. Preference will be given to workshops using DERIVE for WINDOWS as well as 
those using the TI-92. A particular preference will be given to workshops using the CBL or 
CBR in conjunction with the TI-92. Other platforms will be considered depending upon the 
availability of microcomputer laboratory and seminar room space. Attendance at a workshop 
will be based on participants' choices at registration and be limited to 25 in addition to the 
instructor. 

Format for Proposals and Deadline 

Proposals for papers should contain the title of the paper, the names and associations of the 
authors, the name and preferred address of the contact person for the paper, the length of the 
presentation (25 min or 50 min) and an extended abstract of two pages double spaced with 
one standard margin. 

Proposals for workshops should contain the title of the workshop, the name and associations 
of the presenters, the name and preferred address of the contact person, the intended audience 
for the workshop, the level (introductory, advanced or subject oriented) of the workshop, a 
brief abstract, and a detailed outline of the workshop (2 pages total.) 
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Four copies, collated and stapled, should be mailed to 
Professor Carl Leinbach 

Department of Mathematics & Computer Science 
Box 402 

Gettysburg College 
Gettysburg, PA 17325 
USA 

Please mark outside of your envelope: International DERIVE & TI-92 Conference. Papers 
must be received in Gettysburg by November 15, 1997. Notification of acceptance will be 
made during January of 1998. 

Keynote Speakers 

During the morning of three days of the conference there will be an invited keynote address 
for all participants in the conference. No activities will be scheduled during the keynote ad- 
dress. The following individuals have been invited and tentatively accepted our invitation to 
be keynote speakers. 

Josef Bohm Austria 

Wade Ellis USA 

Adrian Oldknow UK 


The Committee will also invite David Stoutemyer of the Soft Warehouse to make a special 
presentation to the conference. 


Conference Committee (Tentative) 

The conference organizers are 
Carl Leinbach 
Gettysburg College 

Department of Mathematics & Computer Science 
Gettysburg, PA 17325 

e-mail: leinbach@cs.gettysburg.edu 


Bert K. Waits 
The Ohio State University 
Mathematics Department 
231 W. 

18th Ave Columbus, OH 43210 
waitsb@math.ohio-state.edu 


all correspondence regarding the Third International DERIVE and TI-92 Conference should 
be addressed to either of above. Proposals are to be submitted to Carl Leinbach by mail or e- 
mail. 


The following people have been invited to serve as members of the Conference Committee 
Barbel Barzel, Germany Terence Etchells,UK 

John Berry, UK Bernhard Kutzler, Austria 

Josef Bohm, Austria Josef Lechner, Austria 
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Marvin Brubaker, USA John Monaghan, UK 

Paul Drijvers, Netherlands Jeanette Palmiter,USA 

Robert Hill, USA Barbara Vincent, USA 

Committee members will assist the organizers with the refereeing of papers and the general 
organization of the conference. 

History of the Conference 

The first International DERIVE Conference was held during the summer of 1994 in 
Plymouth, England and was organized by Professor John Berry of the University of 
Plymouth. The con-ference hosted over 200 participants and set the format for future 
conferences. 

The Second International DERIVE Conference was held during the summer of 1996 at 
Schloss Birlinghoven in Bonn, Germany. This conference was organized by Barbel Barzel. 
This conference hosted more than 250 participants and featured several presentations using 
the TI-92 as well as DERIVE. 

For the Third International DERIVE and TI-92 Conference the TI-92 is given full partnership 
in the title of the Conference. It is hoped that teachers who are using hand held computer al- 
gebra systems will make several contributions to the conference in the form of presentations 
and workshops. 


Registration Fee 

The registration fee for the conference ist expected to be $150.00 US for participants, $100 
for accompanying adults (includes allied program, battlefield tour, and conference trip) and 
$50 for children under 12 (also includes allied program, battlefield tour, and conference trip.) 
This fee will be confirmed when the formal registration forms are distributed. Requests for 
these forms may be made by contacting Carl Leinbach at his e-mail address of 

" leinbach@cs . gettysburg.edu" 

or by mail to him at the address given above for the submission of proposals. 

Housing 

Because Gettysburg is the site of an important battle of the USA Civil War (1860-1865), it is 
a very popular tourist attraction. Enclosed with this call is a listing of hotels/motel accomoda- 
tions in the Gettysburg Area. Those marked with a + are within two blocks walking distance 
of Gettysburg College (the host institution for the Conference). Other motels are an easy drive 
from the College or a 20 - 30 minute walk. The cost of motel rooms varies, but usually runs 
from between $80.00 and $100.00 per day. It should be noted that Gettysburg has several 
restaurants ranging from fast food to fine dining. 
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3D-Grafik mit Adjazenzlisten 

3D-Graphics Using Adjacens Lists 

Hartmut Kfimmel, Biedenkopf, Germany 


1. Einfuhrung 

Zur Darstellung linearer Abbildungen durch Matrizen kann man die Matrix von links oder von rechts 
mit dem abzubildenden Vektor multiplizieren. Durch Transponieren kann eine dieser beiden Schreib- 
weisen in die andere fibersetzt werden - beide sind im Bereich der Computergrafik fiblich. Im Folgen- 
den habe ich mich ffir die zur gewohnlichen Schreibweise bei Gleichungssystemen passende Darstel- 
lung y = A • x entschieden. Dabei kann das Bild einer Figur durch eine einzige Matrixmultiplikation 
berechnet werden, wenn die Punkte als Spalten einer „Figur-Matrix“ aufgefasst werden. 

Die Plot-Funktion von DERIVE erwartet bei Punktlisten die der Reihe nach zu zeichnenden Punkte in 
Zeilen einer „Figur-Matrix“. Diese DERIVE - Syntax ist nicht „kompatibel“ zur oben bevorzugten 
Schreibweise fiir Abbildungen und Gleichungssysteme, bei der die Matrix von links multipliziert wird 
und Punkte als Spalten aufzufassen sind. Zur Darstellung komplexerer Figuren erscheint mir die 
„Punkt-Reihung“ ohnehin nicht so gut geeignet, weil manche Kanten mehrmals „durchlaufen“ werden 
mfissen. 

1 . Introduction 

Linear mappings can be performed by a multiplication of the transformation matrix with the 
vector to be transformed from the left- or right hand side. I chose the form y = A . x - 
according to the presentation of simultaneous linear equations. Then the points are 
interpreted as columns in a "figure-matrix". We have to consider that DE RIVE's PLOT- 
f unction expects the points to be plotted as rows in a matrix. From my point of view this way 
does not seem not to be comfortable for more complex figures, because in most cases some 
edges have to be run through several times. 


2. Reprasentation von Kantenmodellen als Graphen 

Ein Graph besteht aus Knoten und Kanten. Zur Reprasentation eines Graphen ist anzugeben, welche 
Knoten mit einer Kante verbunden sind. Zur geometrischen (grafischen) Darstellung wird zusatzlich 
die Information fiber die Orte der Knoten benotigt. In der Graphentheorie sind hierffir (mindestens) 
zwei Losungen gebrauchlich: Adjazenz-Matrizen und Adjazenz-Listen. 

Das Matrixelement einer Adjazenzmatrix gibt 
an, ob die beiden Knoten mit dem Index i und k 
zu verbinden sind: 

f 1 falls P t mit P k verbunden 
lk [0 sonst 


Furaile Punkte mit Index i=l..n 


Fur alle Punkte mit Index k=l..n 


lstA jk >0? / 
ja / nein 

Zeichne Kante P R 
i k 



Die Darstellung des Graphen erfolgt nach dem 
im Struktogramm angegebenen Verfahren. 
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Eine Adjazenzliste enthalt fur jede Kante ein Paar der Indizes der zu verbindenden Knoten. So eine 
Liste kann beispielsweise durch eine zweispaltige Matrix realisiert sein: 

f A n Index des 1 . Knotens 
Zeile A. =< J ’ 

J I A j 2 Index des 2. Knotens 

Die Darstellung des Graphen erfolgt nach dem 
im Struktogramm angegebenen Verfahren. 

Die Darstellung als Adjazenzmatrix hat Vorteile bei der Suche nach Wegen im Graphen. Bei ungerich- 
teten Graphen steht jedoch zuviel an Information zur Verfugung - die „halbe“ Matrix wiirde geniigen. 
Fur die Darstellung von Kantenmodellen eines Korpers mit den Mitteln von DERIVE erscheinen des- 
halb Adjazenzlisten als geeignete Struktur. 

2. Presentation of edge-models as graphs 

Usually a graph consists of nodes and edges. It must be known which nodes are connected 
by an edge. This information can be expressed by an adjacens matrix or by an adjacens list. 
The element of the matrix A ik = 1, if point P, is connected with P k , otherwise A ik = 0. The list 
consists of index-pairs, which describe all the edges by pairs of numbers, each of them 
belonging to one point. 

3. Beispiel fUr die Darstellung eines Kantenmodells 

Fur ein Pentagon und ein Pentagramm 
ist jeweils die Adjazenzliste angege- 
ben: Alle Adjazenz-Angaben beziehen 
sich nur auf die Indizes der Knoten. 

Zur grafischen Darstellung sind zu- 
satzlich (z.B. in einer Matrix unter 
diesen Indizes) die Koordinaten der 
entsprechenden Punkte anzugeben. 





bzw. in Matrixschreibweise: 


bzw. in Matrixschreibweise: 
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3. An example for presentation of an edge-model and 

4. How to do it in DERIVE 

You can find two lists - one producing a pentagon, the other one will lead us to a pentagram, 
both using the same five points penta, resulting from 

penta := VECTOR([COS(c|>) , SIN(<D) ] ,<D, 0° , 360° , 72°) 

( The point list is given in form of a n x 2 matrix with x- and y-coordinates in two rows.) 

The two lists - matrices with 2 columns - gon and gram are describing how to connect the 
points given in the list of points called penta. The function M0VE( . . . ) helps to translate the 
figure and can be found in the DERIVE-file TOOL3D below. 


4. Umsetzen der Adjazenzlisten-Reprasentation 

Die folgend definierte /)/:/?/ V7:-l'iinkt ion KANTEN erlaubt es, die Kanten einer Figur mit Hilfe einer 
Adjazenzliste zu plotten. 

Dabei wird auch das hier benotigte Transponieren der spaltenweise angegebenen Punktliste in die von 
DERIVE erwartete zeilenweise Auflistung ubemommen: Die Fn nk tion MOVE ist weiter unten angege- 
ben! 



AUTHOR expression: penta : SECTOR C LCOSta », SIMti ) 1, $,0 


Inter expression (press FI for help) 

User G :SALL£HATR\CALG\DF! Free:96>: Lin Define Algebra 


This is the DERIVE screen from 1997. Next page shows the DERIVE 6 screen with 
KANTLIST.MTH and TOOL3D.MTH preloaded as utility files. 
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KANTLIST.MTH 

#1: [Precis ionDig its := 12, Precision := Approximate] 

#2: [NotationDigits := 4, Notation := Decimal] 

#3: KANTEN-Tool zur Darstellung von Figuren mit Adjazenzlisten 

#4: Eckpunkte stehen in SPALTEN von fig, M 1 ist die Transponierte von M 

#5: Vertices are given in columns of fig, M 1 is transposed of M 

KANTEff lg , ad j , i ) := |"f ig 1 , fig 1 "I 

#6: L adjiiil ad j a i 1 2 J 

#7: KANTENff ig , adj) := VECTOR(KANTE(f ig , adj, i), i, 1, DIM(adj)) 

#S: Nun Figuren (Punkte spaltenweise) und Adjazenzlisten eingeben 


*** Enter figure (points in columns) and ad jacens-1 ist *** 


5. Zur Verstandlichkeit des KANTEN-Tools im MU 

Die Funktion KANTE verwendet das Konzept der indirekten Adressierung. Das ist eine elegante, fur die 
Informatik typische Denkweise, zugleich aber auch eine Anwendung von Matrizen. Wie aufwandig ist 
es, die Arbeitsweise des KANTEN-Tools durchschaubar zu machen? Im Unterricht wiirde ich die 
Auswertung des Ergebnis an einem Beispiel-Aufruf durchspielen: 


Bestimme die durch KANTEN(quad , umf) be- 
rechnete Figur (#8): 

#8: 


quad 


' 7 10 

. 1 1 


10 7 

4 4 J 


umf := 


1 2 
2 4 

2 3 

3 4 
L 4 1 
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Die schrittweise Auswertung von KANTE(quad , umf , 2) ergibt: 


#9: [quad :=, umf :=] 

KANTEfquad, umf, 2) = 


# 10 : 


j^quad 1 


mf j, 2 j,l 


, quad 1 


#11: umf := 


#12: 


1 2 
2 4 

2 3 

3 4 

4 1 

KANTEfquad, umf, 2) = 


j'quad 1 , quad 1 j 


umf i2i2_ 


#13: 


#14: 


quad := 


7 10 10 7 
L 1 1 4 4 J 


KANTEfquad, umf, 2) = 


10 

7 


1 

4 


1. Auswertung : Ausdruck #12 (quad noch nicht definiert!) 

2. Auswertung: Ausdruck #14 

Nun folgt die Konstruktion einer Liste solcher Punktepaare: 


#13: [quad :=, umf :=] 

#16: KANTENfquad, umf) = VECTORf KANTEfquad , umf, i), 


#17 : umf 


1 2 ' 
2 4 

2 3 

3 4 

4 1 . 


#18: KANTENfquad, umf) 


quad 1 

quad 1 

1 

2 

quad 1 

quad 1 

2 

4 

quad 1 

quad 1 

2 

3 

quad 1 

quad 1 

3 

4 

quad 1 

quad 1 


4 1 J 


1, DIMfumf)) 
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#19: quad 


' 7 10 10 7 
. 1 1 4 4 


#20: KANTENCquad, umf) 



1 

1 


10 1 
7 4 


10 1 
10 4 


' 10 4 

. 7 4 


7 

7 


4 

1 


I 


Hier soil nur die Arbeitsweise einer vorgelegten Funktion nachvollzogen - und verstanden - werden. 
Die Fahigkeit, solche Funktionen in weniger vertrautem Zusammenhang zu programmieren, ist ein 
vollig anderes Lemziel! 


5. How to explain the tool in Math Teaching 

The function kante uses the concept of indirect addressing. This is an elegant way of 
thinking which is typical in information technology, and it is an application for matrices. 

I believe that it might make sense to make this tool understandable even in standard 
mathematics teaching. Expression #8 from above describes a special figure. How does it 
look like? 

I would demonstrate the concept for stepwise evaluating one call of kante. In this case I want 
the students to understand how a ready made function is working. Developing such functions 
is quite a different learning goal. 


6. Raumliche Abbildungen 

Als Projektion zur Darstellung von Schragbildem wurde die senkrechte Projektion auf die 1-2-Ebene 
festgelegt, um nicht in der Vielfalt der Moglichkeiten unterzugehen. Die Matrizen fur diese Projektion 
und fur die Drehungen um die Koordinatenachsen sind Standardbegriffe. Zum Erzeugen der ge- 
wiinschten Ansichten muss die Urfigur passend abgebildet werden! 

Die Verschiebung wird hier als spezielle, nicht-lineare Abbildung eingefuhrt - etwas umstandlich 
wegen der Darstellung der Punkte in Spalten der Figur-Matrix. Man konnte sie auch unter Verwen- 
dung homogener Koordinaten in die Standarddenkweisen „einbetten“. 


#1: TG0L3D: Hilfsmittel fur 3D-Grafik 

#2: [Precision := Approximate, Precis ionDigits := 12] 

#3: [NotationDigits := 4, Notation := Decimal] 




'10 O' 

#4: 

DREHIO) 

0 COSO) - SINO) 



. 0 SINO) COSO) . 



' COSO) - SINO) o ' 

#5: 

DREHBO) 

SINO) COSO) o 



0 0 1. 


MOVEff, v) := VECTOR (VECTOR (f + v , 

z, i z 



' COSO) 

0 

SINO) ' 


DREHZO) 

0 

1 

0 



. - SINO) 

0 

cosO) . 



proj!2 


10 0 ' 
. 010 . 


i, 1, DIM(f )), z, 1, DIMff)) 
1 


#6: 
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Die Darstellung des Einheitswiirfels zeigt den Ubergang zur dritten Dimension. 

Das KANTEN-Tool bleibt unverandert, da die Adjazenz-Information auch fur die projizierten Figu- 
ren gilt! 


#1 : L0AD(H : \DOKUS\DNLs \DNL97\WTH26\KANTLIST . MTH ) 

#2 : L0AD(H : \D0KUS\DNLs\DNL97\MTH26\T00L3D. WTH) 


0 110 0 1101 


#3: unit_c 


00110011 


L0 0 0 0 llllj 


#4: c_adj 


123456781234' 

234167855678. 


6. Mappings of 3D objects 

In order to produce oblique views I chose the normal projection onto the 1-2 plane. The 
matrices for this projection and for the rotations about the coordinate axes are standard. 
Translation is introduced as a special non-linear mapping. You might prefer using 
homogeneous coordinates. 

The presentation of the unit cube shows the change to the 3rd dimension. The KANTEN-tool 
will remain the same because the adjacens list is also valid for projected figures. 



7. Ein Abschluss-„Film“ ... 

Die Herstellung des Films als Liste von Einzelbildem zeigt die typische, daten-orientierte („objekto- 
riente“) Arbeitsweise im funktionalen Stil: 

Als Grundobjekt wird ein „Turm“ durch seine Punktliste tower und die Adjazenzliste tadj definiert. 
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Durch Abbildungen wird zuerst jedes „i-te“ Einzelbild und dann daraus der Film als Liste der Einzel- 
bilder hergestellt. 

7. Finally a "Movie".... 

The "tower" is defined by its point list tower and its adjacens list tad j . The pictures are the 
result of a rotation combined with a translation. The movie consists of a list of single pictures. 

(DREHl(ot), ... are identical to ROTATE_X(oO,... from GRAPHICS. MTH.) 


# 1 : 

# 2 : 


#3: 


#4: 


LGAD(H: \DOKUS\DNLs\DNL97\WTH26\KANTLIST. WTH) 

LOAD (H:\DO KU 5 \ D N L s \ DN L 97 \ MTH 2 6 \T00 L 3 D . MTH ) 
"01100110 0.5 

00110011 2 
.00001111 0.5 

123456781234347 
234167855678999 


tower := 


tad j := 


8 

9 


#5: tow(ffl) := DREHl(cx) ■ DREH2(30 ■ 1°) ■ tower 

#6: tower_pic(i) := M0VE(tow(15 ■ i ■ 1°) , [0.75-i, 2, 0]) 


VECTORCKANTENCproil2.tower_picC2.i), tadj), i, 0, 7) 



In DERIVE 6 you can really have a movie using a 
slider bar! 

KANTEN(projl2 ■ tower_pic(2 ■ i_) , tad j) 



Introduce a slider for i_ with 0 < i_ < 7 and 70 
intervalls. 

Then drag the slider and observe the tower moving 
smoothly across the screen. 



Josef 
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8. Adjazenzlisten fiir orientierte Flachen 


AbschlieBend mochte ich zeigen, dass der von Hubert Weller im nachsten DNL erklarte Sichtbarkeits- 
test mit orientierten Flachen „nahtlos“ in das Adjazenzlistenkonzept eingebunden werden kann. 

Zur Darstellung eines Korpers miissen zunachst die Punkte spaltenweise eingegeben werden. Als Ad- 
janzenzlisten-Eintrage fiir jede Flache sind die Indizes der beteiligten Punkte so einzugeben, dass die 
Reihenfolge einen Linksumlauf beschreibt, wenn die Blickrichtung auf das Zentrum des Korpers 
zugrunde gelegt wird. Aus diesen Daten erzeugt die Funktion f 1 as eine Liste von Flachenlisten zur 
Darstellung durch die Plot-Funktion, aus der die sichtbaren Objekte mit show und die unsichtbaren mit 
unshow ausgewahlt werden konnen. 


8. Adjacens lists for oriented surfaces 

At the end of this article I’d like to demonstrate that Hubert Weller’s visibility test (see next 
DNL) can be implemented into the concept of adjacens lists. 

For presenting a figure one has to enter the points as columns. The entries for the adjacens 
lists of the single surfaces are the numbers of the points involved together with the sequence 
of the indeces of the points in an order describing a left hand circulation when viewing at the 
center of the body. Function f 1 as creates a list for plotting, show is collecting the visible 
objects while unshow is collecting the invisible ones. 

Utility file f 1 adl i st . mth 


#1: Verdeckte Kanten - Hidden Lines 

#2: [Precision := Approximate, Notation := Decimal] 

#3: 1. Matrices for linear mappings in R3 


#4: 


id := 


10 0 
0 10 
L 0 0 1 


#3: Translation of object f by vector v 

movetf, v) := VECTOR (VECTOR (f + v , 

#6 : z , i z 

#7: Similar to isometric projection 

#S : iso := R0TATE_X(30 o ) ■ ROTATE_Y(35 G ) 

#9: Vertical projection = top view 


#10: projlZ 


10 0 " 
0 10 . 


DIMff )), z, DIMff )) 
1 


#11: 2. Special vector operations: project and transpose 

pr!2(lst) := VECTOR (projl2 -1st , k, DIM(lst)) 

k 


#12: 
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transp(lst) := VECTORflst ’ P i P DIMflst)) 

#13: i 

#14: 3. Producing a list of areas from points and adjacens data 

flafpts, 1st) := VECTOR(pts 1 P i P DIM(lst)) 1 
#15: lst|i 

flas(pts P fad) := VECTOR(f laCpts p fad ) P j p DIMffad)) 

#16: j 

#17: 4. Visibility check of oriented surfaces using determinants 

showdet(lst) := DET([lst - 1st p 1st - 1st 1) 

#18: L 1 23 2j 

show(flist) := VECTQR(IF(showdet(f 1 ist ) < Q P flist ) P j P DIM(flist)) 
#19: j j 

unshowff 1 ist) := VECT0R(IF(showdet(f 1 ist ) > 0 P flist ) P j P DIM(flist)) 

#20: j 3 

Nun sind die Korper iiber Punkte und orientierte Flachen zu definieren. 

Now the bodies can be defined by points and oriented surfaces. 


9. Orientierte Flachen und Sichtbarkeit bei konvexen Korpern 

Der folgende Film wird aus einem „Turm" mit passenden Flachen- Adjazenzlisten erzeugt. Zur De- 
monstration ist die Definition des Korpers angegeben. 

Die kompletten DERIVE-Arbeitsblatter finden sich auf der Homepage. Teilweise miissen Aufrufe von 
Funktion vor dem Plotten ausgewertet (vereinfacht bzw. approximiert) werden, um die Daten fur die 
Plot-Funktion zu erzeugen! So wird Platz fur umfangreiche Punktlisten eingespart. 

9 L Oriented surfaces and hidden lines of convex bodies 

The following “movie” is created from a “Tower” with appropriate surface-adjacens lists. 

You will find the complete DERIVE-listings on the home page. Some functions must be 
simplified before plotting to create the plot data. 

#1: L0AD(H : \D0KUS\DNL\DNL26\f ladl i st . mth) 

#2: Once more the tower 


[01100110 0.5 

tower := 00110011 2 

Lo 0 0 0 1 1 1 1 0.5 


#4: fadj := [[1, 4, 3, 2, 1], [2, 3, 7, 6, 2], [1, 2, 6, 5, 1], [1, 5, 8, 


4, 1], [5, 6, 7, 8, 5], [3, 4, 9, 3], [3, 9, 7, 3], [7, 9, 8, 7], 
[4, 8, 9, 4]] 
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#5: 

#6: 

#7: 

#3: 

#9: 

#10: 

# 11 : 

#12: 

#13: 


Generating the surfaces for the i-th picture of the movie 
tow(oc) := ROT ATE_X (w ) ■ ROT ATE_Y ( 3 0 ° ) ■ t owe r 

3- 


towerpic(i) := move tow(15 ■ l °) , 


f 2 ' ° 1 ) 


towersurf(i) := transp(prl2(f las (towerpic(i ) , fadj))) 
Test: Plot the next expressions! 

test_tower := transp(prl2(f las(towerpic(2) , fadj))) 
view_ tower := s hi ow ( t e s t_t owe r ) 


VECTOR(show(towersurf (2 ■ j)) p j, 0, 6) 



Plot of #10 


Plot of #11 



The Movie 
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As you might know from other DNLs, I am very keen in working with graphic representations. 
So it was a "must" for me to use Hartmut s Toolbox. 

I tried to produce my own movie with a figure having a hexagon as its base and a square as 
its top. 


#46: 


#47: 


#48 : 
#49: 
#50: 


body := 


badj := 


4 2 -2 -4 -2 2 2 -2 

0 2-J3 2-J3 0 - 2-J3 - 2-J3 2 2 

00 00 0 055 

12345678 9 10 123 

23456189 10 7778 

BOD(h) := ROTATE _X ( « ) ■ R0TATE_Y f 3 0 ■ 1 " ) ■ b od y 
BODPIC(i) := M0VE£B0D(15 ■ i ■ 1°) , [5-i, 4, 0]) 
VECTOR(KANTEN(pro jl2 ■ B0DPIC(2 ■ i ) , badj), i, 0, 7) 


-2 2 ‘ 

-2 -2 
5 5 . 

4 4 5 6 

8 9 9 10 


1 ' 
10 . 



#51: fb := [[1, 6, 5, 4, 3, 2, 1], [7, 8, 9, 10, 7], [1, 10, 6, 1], [1, 7, 10, 

1], [2, 7, 1, 2], [2, 3, 8, 7, 2], [4, 8, 3, 4], [4, 9, 8, 4], [4, 5, 9, 
4], [5, 6, 10, 9, 5]] 


#52: bsurf(i) := transp(prl2(flas(B0DPIC(i) , fb))) 

#53: VECT0R(show(bsurf(2. j)), j, 0, 7) 



In the next DNL you will find Hubert Weller's contribution dealing with the hidden lines. This 
"geometric" series" will be completed by some of my own tools for special projections 
(axonometric and central projection, ....) 

It might be attractive to compare the various approaches. Josef 
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DERIVE for DOS and DERIVE for WINDOWS 

Version 4.05 


Dear DERIVErs, 

On 5/12/97 Soft Warehouse, Inc. released DfD and DfW 4.05. New functions have been added to 

existing utility files and more contributions have been made to the user directory. 

DfW 4.05 makes the following interface improvements over 4.04c: 

1. Definite integrals, sums, and products (e.g. INT(SIN(x)/x, x, 0, t) ) can now be plotted directly in 
the 2D-plot window without simplifying or approximating first. 

2. When in Case Sensitive Input Mode both upper & lowercase variables can be selected in Calculus 
and Solve Algebraically listboxes. 

3. Horizontal scrolling has been adjusted so that the smallest scrolling increment is one character. 
Additionally, horizontal scrolling is no longer dependent on the selected expression. The Home and 
End keys now scroll all the way to the left and right, respectively. 

4. Parametric plots now use the given minimum and maximum parameters (in degrees) when in 
’’Degree” angle mode. 

5. If-statements using the Boolean operator NOT can now be plotted. 

6. Allow vectors of the form [z=2-x-y] to be written to an AcroSpin file instead of having to select the 
equation within the brackets. 

7. The Online Help has been updated; additional indeces and topics have been added. New topics 
include: Animation with AcroSpin ’97, Using and Saving Algebra State Variables; DERIVE 
Initialization File; and Windows 3. lx Installation Questions. 


Soft Warehouse, Inc. (the authors of DERIVE. A Mathematical Assistant) 

3660 Waialae Avenue, Suite 304 

Honolulu. HI 96816-3259 U.S.A. Tel: (80S) 734-5801 Fax: (808) 735-1105 

Emai 1 : swh@. aloha . c om Web page: http : // w w w . derive . c om 


Shading regions - 
- See page 43!! 
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Herr Peter Graser hat den Wunsch nnit einer osterreichischen Schule in Kontakt zu treten, unn even- 
tuell gemeinsame DERIVE - Projekte durchzufuhren. Gerne drucke ich seinen Brief hier ab und wurde 
mich freuen, wenn uber den DNL eine Zusammenarbeit zustande kame. Josef 

ich arbeite am Willi-Graf-Gymnasium Saarbriicken. Saarbriicken ist die Hauptstadt des Saarlan- 

des, eines kleinen Bundeslandes von Deutschland an der franzosischen Grenze. 

Vor einigen Jahren fuhr ich im Auftrag der Fachkonferenz Mathematik/Physik/Informatik nach Wien 
um mich uber den Einsatz von DERIVE im Unterricht kundig zu machen. Dort informierten mich 
Herr Dr. Heugl und Herr Klinger aufs ausfuhrlichste. 

Da unsere Schule bald einen NetzanschluB haben wird, suche ich eine Schule in Osterreich, die an 
gemeinsamen Projekten (vor allem in Mathematik mit DERIVE - Einsatz) uber das Netz interessiert 
ist. Bitte helfen Sie mir eine solche Schule zu finden. 

Mit freundlichen GriiBen, 

P. Graser p.graeser@hit.handshake.de 


Die zweite Information betrifft ein neues Diskussionsforum: 

From: Erich Neuwirth, neuwirth@smc.univie.ac.at 
Liebe Kollegen! 

Da in Osterreich Computer- Algebra ziemlich weit verbreitet im Mathematikunterricht eingesetzt wird, 
scheint es mir vernunftig zu sein, ein diesem Thema gewidmetes Diskussionsforum zu schaffen. daher 
gibts seit heute auf dem listserver auf sunsite.univie.ac.at eine neue mailing list 

cau-1 (computer-algebra im unterricht) 

Das Fomm soil dem Austausch von Ideen zu diesem Thema dienen und unter anderem erfolgreiche, 
aber auch weniger erfolgreiche Projekte beschreiben. Beziiglich der eingesetzten Systeme gibts keine 
Einschrankungen, Derive, Mathematica, Maple .... sind alle gleich willkommen. 

sunsite austria hat die Sachwidmung: Computerunterstiitzung in mathematischer, statistischer und 
naturwissenschaftlicher Ausrichtung. Daher werde ich auch versuchen, auf der sunsite Sammlungen 
von Beispielen und Programmen fur Computer- Algebra- Systeme anzulegen. Die neue mailing list soil 
daher auch dazu dienen, Informationen, wo auf dem Internet Resourcen uber Computer- Algebra ver- 
fiigbar sind zu sammeln. 

Wie man auf der neuen Liste computer-algebra im unterricht subskribiert: man schickt mail an 
listserv@sunsite.univie.ac.at 
mit dem Inhalt (/nicht/dem Subject oder Betreff) 
subscribe cau-1 <Ihr Name> 

wobei natiirlich <Ihr Name> durch Ihren Namen ersetzt werden muss. Die Liste kann von alien Leuten 
subskribiert werden. 

Bitte subskribieren Sie zahlreich, 

Erich Neuwirth 

PS.: Normalsprache auf dieser Liste sollte deutsch sein. Allerdings soil das niemanden davon abhalten, 
Beitrage, die auf dem Netz gefunden wurden, aber beispielsweise englisch sind, auch uber die Liste zu 
verbreiten. 

Visit my Spreadsheet home page at http : / /sunsite . uni vie . ac . at /Spreads it e 
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Comment of 201 1: Some of the functions created in this article by "Mr. Titbits”, Johann 
Wiesenbauer, have been implemented in DERIVE since 1997. Hence all of Johann’s 
complaints are of “historical” interest. It must be mentioned that Johann contributed 
enormously for improving the number theory functions of DERIVE through many years, Josef 


Titbits from Algebra and Number Theory (10) 

by Johann Wiesenbauer, Vienna 

This time Tm going to deal with a topic that I have put off several times because it scares 
me a little bit. It is the present state of the utility file number.mth. What gives me the jit- 
ters is the question whether I should adopt an attitude that everything is in apple-pie or- 
der (apart from, you know, minor insufficiencies that hardly deserve to be mentioned!) or 
risk the highly esteemed friendship with the people at the Soft Warehouse by not mincing 
matters. Well, as you may conclude from the last sentence I have made up my mind to 
choose the latter option and this is what I think of the utility file number. mth (referring to 
the current version 4.04 of DERIVE). It s a real mess consisting of functions that range 
from being very useful to totally useless, or to put it in another way, from being terrific to 
terrible. 

"Terrible? What the hell is he talking about?’ 1 , you might be asking. Well, let me give you 
some examples. At first two extremely simple functions where the function call is longer (!) 
than the definition of the function: 

SQUARE J/VAVE(x):=(-1 ) A FLOOR(x) 

POWER_MOD(n l d J m):=MOD(n A d l m) 

Okay, okay, Tm petty and overcriticah I know. I promise to refrain from further examples 
of this sort. But what about the following function INVERS E_MOD ( m 5 n) that is supposed 
to return the inverse of m modulo n? 

INVERSE_MOD(m l n):=IF{m=1,1 ,IF(m=n 5 ?, 

(n*INVERSE_MQD{m-MOD(n.m),m)+1)/m)) 

As a first test I wanted to calculate the inverse of - 1 mod 10007 . But after simplifying 

INVERSE_MOD(-1 ,10000) 

the result was a question-mark (!). Well, obviously they didn’t take into account negative 
numbers, but this is only a small inconvinience as it can be easily circumvented by simpli- 
fying the equivalent expression 

INVERSE_MOD(9999,10000) 

This time I got the message "Insufficient computation memory” after some seconds. (By 
the way, my Pentium 200 PC has got 32 MB RAM!) After some more attempts I found out 
that this routine only w-orks if the modulus n doesn't exceed a number in the range of sev- 
eral thousands. Very strange, indeed! After all. I had shown them a number of times (in 
DXL #14. #17, #24) how it is done properly! Those routines of mine based on Euclid's Ex- 
tended Algorithm (EEA) will yield the correct result within fractions of a second even if n 
has 100 digits or more (Check it!). 

I hate to say it. but the next routine SOLVE_MQD(u,x,n) which is supposed to solve linear 
conguences of the form u(x) = 0 mod n is an even bigger disaster, both in terms of pro- 
gramming and w T hen it comes to performance. 
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SOLVE_MOD_AUX3(s J d ? m):=VECTOR(s+m*k_ J k_ l O,d-1 } 

SOLVE_MOD_AUX2(c.d l m J x J n):=IF(MOD(c ) d)=Q,SOLVE_MOD_AUX3( 

MOD{c/d*INVERSE_MOD(m/d P n/d} J n/d} ! d,n/d) J ?) 

SOLVE_MOD_AUX1{u ! m T x J n):=SOLVE_MOD_AUX2(m*X"U l GCD(m 1 n) J m P x l n) 

SOLVE_MOD(u f x,n):=SOLVE_MOD_AUX1(u l DIF(u r x),x J n) 

As you can see it uses 4 (!) auxiliary functions including the notorious INVERSE_MGD. In 
principle, my attitude towards auxiliary functions is the following: Avoid them unless they 
are unavoidable (as it is often the case with recursive programming) or unless you want to 
make the program more readable for others, i.e. for didactic reasons. Sometimes an auxil- 
iary function might also be interesting on its own (IXVERSE_MOD would be of this kind 
if it worked properly). In all other cases the abundant use of auxiliary functions (like 
above!) is usually a sign that the programmer was either not able or too lazy to create a 
self-contained program. 

As for performance it comes as no surprise that SOLVE_MOD inherits all the weaknesses 
of IXVERSE_MOD. which means that it is rather slow and works only for small n. On top 
of it all. unsolvable congruences will yield a question-mark instead of the empty list [ ] ! 

How to fix this problem? Well, in the first place we need a working IXVERSE_MOD. I 
would have liked to avoid it. believe me, but here it is (for the fourth time!): 

INVERSE_MOD(a,m):=IF(GCD(a,m)=1 5 MOD((ITERATE(IF(MOD(a_ l bJ=0,[a_.b_ l c_ 5 dJ 
, [b_ 1 MOD(a_,b_) : d_,c_-FLOOR(a_ : b_)*dJ) : [a_ 1 b_ J c_.d ,0])) SUB 4jm)) 

And this is - for purists only! - my version of SOLVE_MOD without any auxiliary functions 
apart from INVERSE_MOD (u is supposed to be either a linear term or a linear equation 
in the variable x): 

SOLVE_MOD(u,x,m}:=ITERATE(IF(MOD(LIM(RHS{u)-LHS(u) p x 1 0) 5 t_)=O, 

ITERATES(s_+m/t_,s_, MOD(INVERSE_MOD(DIF(LHS(u) J x)/t_ 1 m/tJ* 
LIM(RHS(u)-LHS(u) s x J 0)/t_ l m/t_) f t_-1) J [ ]) P t_ l GCD(DIF(LHS(u) 1 x) J m} T 1) 

It goes without saying that this routine won't let you down even when it comes to numbers 
with several hundreds of digits and more (Check it!) 

One more example concerning the function PELL(n) which returns the n-th number in the 
Pell sequence defined recursively by 

1.^=0 

Some time ago (cf. my 'Titbits' in DNL.#21) I showed them how to program the Fibonacci 
sequence properly by means of the Lucas sequence. That program was polished up by A1 
Rich himself (cf. DXL #22. p.33) and was used to turn the former "terrible” version of 
FIBGNACCI(n) in number. mth into a "terrific" one. 

Well, you won't be surprised, if I tell you that a similar trick also works in the case of the 
Pell sequence. All you have to do is to write a program that computes the numbers in the 
generalized Lucas sequence defined recursively by 

> L 0= 2 i L l=P 

where the integers p and q are the parameters of the sequence. Just to be on the safe side 
we ca rry out these calculations mod m (note that m=0 leaves the numbers unchanged, if 
it is that what you want!) which leads to the following programming code 
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GEN_LUCAS(n ) p.q.m):=IF(n=0 1 MOD(2,m),(ITERATE(IF((n AND d_)=Q, 
[MOD(a_ A 2-2*MOD(q A c_ l m) ! m) f MOD(a_*b_-p*MOD(q A c_ J m) l m) J 2*c_ T d_/2] J 
[MOD(a_*b_-p*MOD(q A c_ J m) T m) J MOD(b_ A 2-2*MOD(q A (c_+1 ) J m} l m) J 
2*c_+1,d_/2]} J [a_,b_ J c_ l dJ J [2 J p f O l 2 A FLOOR(LOG{n p 2))] 1 FLOOR(LOG(n l 2))+1))SUB 1) 

LUCAS{n):=GEN_LUCAS(nJ ,-1 ,D) 

FIBONACCI(n):=(2*GEN_LUCAS(n+1 ,1,-1 r 0)-GEN_LUCAS(n,L-1 ,0}}/5 

PELL(n):=(3*GEN_LUCAS(n,2 f -1 ,0)»GEN_LUCAS(n+1 ,2,-1 ,Q»/4 

PELL(1 0000)=836..665 (3827 digits; 0.1 sec vs. 4.5 sec with the "old” routine!) 

Just to give an example where in is not 0, here is a nice probabilistic primality test based 
on generalized Lucas sequences: 

LUCAS_PRIME(n ? p}:=IF(GENJJJCAS(n,p,1,n)=p, true, false) 

LUCAS_PRIME(1 0 A 1 00+267, 1 )=true (0.3 sec !) 

I leave it to you to experiment a little bit and find composite numbers that are "pseudo- 
prime" with respect to this test for some choice of p. 

Would you believe that there are functions in number. mth that are never used in the lit- 
eral sense of the word because they are overwritten during the loading process? There are 
even 3 (!) examples of this kind, namely the first definitions of MERSENXE(n). PARTS(n) 
and DISTIXCT_PAKTS(n). (In the case of the partition functions this is all the more inter- 
esting. as their second definitions are incredibly powerful leaving even the built-in func- 
tions of Mathematic a and MAPLE in the dust.) 

Well. I could go on and on pointing out inadequacies, but actually this would miss the 
point. As a matter of fact, the real weaknesses of number. mth are not the functions which 
are included in it. but the functions that are supposed to be there but aren’t (like Euler's 9 - 
f unction. Mbbius' p- function, the divisor functions a and l Jacobi- symbols (a/n) and what 
have you. Needless to say that all these functions and many more have been defined in 
this column!). 

Well, while writing these lines version 4.05 of DfW has just arrived. Will there be big 
changes in the untility-file number. mth making this article obsolete? After all. Josef is 
urging me to deliver it... No. all the blunders described above are still there, but a lot of 
new functions have been added that nobody has asked for. The input of FIBONACCI(O) or 
LUCAS(O) still leads to a very nasty lock-up! Are there any other users of number. mth 
apart from me. If so. what do they think of it? I'll fry to find the answer to this question at 
the next DERIVE- conference in Sweden. See ya! (j.wiesenbauer@tuwien.ac.at) 


Revising the next contribution was a lot of work. The DERIVE-functions from 1997 (Derive for 
DOS, version 3) did not work because of several reasons. One of them is that the SOLVE 
command was accomplished by SOLUTIONS in later versions. Another reason is that the 
first function REL_AUX2(w) does not work and then all further functions which depend on it 


What to do? Rewriting them and using DERIVE 6’s programming facilities. 


I reprint the old version and then append the DERIVE 6 file. 
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From Linear Inequalities to 
Linear Programming (1) 

Josef Bohm, Wiirmla, Austria 

I have been concerned with Linear Programming for a long time. Some time ago I produced 
a DERIVE utility for applying the SIMPLEX-method. I'll refer to this tool at the end of this se- 
quence of contributions (International DERIVE Journal, voll, number 3). When you introduce 
Linear Programming you will start with problems which can be solved graphically. So I tried 
to use DERIVE for that approach. The graphic representation of linear inequalities is one of 
the im-portant basic skills connected with graphic solutions of optimisation problems. 
Unfortunately I didn't know the "griddy" trick with implicit functions - 1 could use the poor 
excuse that I pro-duced this tool at times of version 2.x - so I developed functions to shade 
half planes according to the given inequalities. 

I used Sergey Biryukow's file from DNL#22 to handle the inequations. Unfortunately I am not 
so skilled in programming powerful tools without using auxiliary functions like Johann Wie- 
senbauer, so I apologize for my "programming style". 

#1: [InputMode :z Word, flotation :z Decimal, NotationDigits 3] 

#2: R El AUX2(w) :i 100-IF(w - (LHS(w) = RHS(w))) + 10-IF(w - (LHS(w) + 1 = RHS(w))) + 

IF(w - (LHS(w) - 1 = PHS(w))) 


' ( 100 = 

1 > 

10 < 

R El AUXl(x) :i SELECT x = rel_aux_ , rel_aux_, 

#3: 1 101 >= 

110 <= 

, l 11 /= 

1,2 

REL(w) :z 

If NUMBER (R El AUX2(w)) 

#4: R El AUX1(REI AUX2(w)) 

■■ ■■ 


SGN(z) := 

If Z > 0 
1 

#5: If z < 0 

-1 
1 

#6: V(u) := DI M ( VAR I AB L ES ( u ) ) 

HELP(u) := 

If V(u) = 2 v (VARIABLES(u)) jl = y 
1 
0 
0 



#7: 
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#8: [TEST(u) := SOLVE(u, y), TEST2(u) := SOLVE(u, x)] 

UIIH(u) :z 

If HELP(u) = 1 
#9: TEST(u) 

TEST2(u) 

AUX(u) := 

If REL(UIIH(u)) = "<z" v REL(UIIH(u)) z "<" 

#10: -1 

1 
1 

d 

#11: S L ( u ) := — RHS(UNH(u)) 

dx 

#12: G(u) := S01UTI0NS( LHS(u) z RHS(u), y) 

P(u) := rx_, 1 i in (G(u)) 1 
#13: L x-*x_ lj 

#14: X_C(u) :z SOLUTIONS(u, x) 

LIM(u) :z 

If HELP(u) z 1 
#15: y = (G(u))U 

x z X_C(u) 

#16: SHV(u , a, b, inc) :z APPEND([LIN(u)] , VECTOR([P(u), [x_, - 100 -AUX(u)] ] , x_, a, b, 

inc)) 

SHH(u , a, b, inc) :z APPEND( [LIN(u)] , VECTOR(rP(u), riOO-SGN(SL(u))-AUX(u), 

#17 : L L 

( p (iO)2 J J • x -> a - b - inc)) 

SW(u, a) :z 

If Sl(u) < 0 a (« < 0 a SL(u) < TAN ( « • 1° ) ) 

-1 

#18: If SL(u) > 0 a (of > 0 a SL(u) > TAN(of-l 0 )) 

-1 

1 

1 


#19: SH(u , a, b, inc, c<) :z APPEND [LIN(u)] , VECTOR IF HEEP(u) z 1 , fP(u) , [ x_ + SW(u, 



SHV(u,a,b,inc) vertical shading for u with a < x < b with increment inc in x-direction 
SHH(u,a,b,inc) horizontal shading for u with a < x < b with increment inc in y-direction 

SH(u,a,b,inc,a) shading for u with a < x < b with increment inc in x-direction showing an angle 
of -tt/2 < a < tt/2 (with +x-direction) 
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Problem for pupils: 

Describe the shaded area by a 
system of inequalities and check 
your proposal. 


This was the original contribution from 1997. 

Yes, I know, that we can now shade parts of the xy-plane described by inequalties. 

x - y > 2 a x + 3<y - 6 < 0 a y > -2 



-2 






-2 



Jm 



10 

12 

14 


* _ 







But it is my strong opinion that in many cases it would be useful demonstrating how this 
region (here the triangle) is created step by step adding one inequality after the other. Doing 
this with DERIVE 6, each plot of an inequality covers the previous ones. 


How can we do it in DERIVE 6? 


REL(w) := 

If (STRING 00)4,1 = v (STRINGOO)j.l = 
# 1 : 1 

-1 

g(u , u_) := 

Prog 

u_ := SOLUTIONS (LHS(u) = RHS(u), y) 

#2: If u_ = [] 

[SOLVE(LHS(u) = RHS(u), x), ®] 
[u_ 4 ,l, d(u_ 4 ,l, x)] 

p(u, x_) := fx_, 1 -im (g(u, u_)) 1 
#3: L x -+x_ lj 
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shh(u, a, b , inc, sh) := 

Prog 

If CgCu, u_))j,2 > 0 v CgCu, u i_))j,2 < 0 

sh := VECTGR([p(u, x_) , ECpCu, x_))il - REL(u) -1000 , Cp(u, x_))-l2]], x_, a, b, inc) 

If CgCu, u _))|2 = 0 
Prog 

If REL(u) = 1 

#4: [a := MIN (a, b) , b := CgCu, u_))j,l] 

[a := CgCu, u_))il, b := MAX(a, b)] 

DISPLAYCEa, b]) 

sh := VECTOR (y = CgCu, u_))|l - i ■ RELC u) - inc, i, 0, (b - a) /inc) 

If (g(u, u_) ) j, 2 = cs 

sh := VECTOR CE RHSCCgCu, uj)|l), x_; RHSCCgCu, u_))-lD - RELCu)-1000, x_] , x_, a, b, inc) 
sh -APPENDCECgCu, u_))ilL sh) 

shv(u, a, b, inc, sh) := 

Prog 

If CgCu, u_))i2 > 0 v Cg(u, u_))i2 < 0 

sh :: VECT0R([p(u, x_), ECpCu, x_))il, CpCu, x_))i2 - SIGNCCgCu, u_))i2) -REL(u) -1000]] , x_, a, 
If (g(u, u_))i2 = «. 

Prog 

If REL(u) = 1 

#5: [a := MIN(a, b), b := RHSCCgCu, u_))il)] 

[a := CRHSCgCu, u_)))ll, b :z MAXCa, b)] 

DISPLAYCEa, b]) 

sh := VECTORCx = a + l-inc, i, 0, Cb - a)/inc) 

If CgCu, u_))i2 = 0 

sh := VECTOR CEx_, CgCu, u_))il; x_, CgCu, u_))il - RELCu)*1000], x_, a, b, inc) 
sh := APPENDCECgCu, u_))il], sh) 


b, inc) 


#7: shhfx - y > 2, -12, 12, 0.25) 

#S: shvO + 3«y - 6 < 0, -4, 13, 0.4) 

#9: sho(y > -2, -S, 14, 0.4, 35) 
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You can distinguish the regions described by the three inequalities. 
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#10: x - y > 2 a )( + ly - 6 < 0 a y > -2 



I superimpose the Boolean expression combining the given inequalities. 

Here are some more examples to demonstrate my function sho: 
sho (x <2, -4, 4, 0.25, 30) sho(x <2, -4, 4, 0.25, 30, -1) 




sho (x + 2-y < 0, -4, 4, 0.25, 45) sho(x + 2-y < 0, -4, 4, 0.25, 45, -1) 




You can enter the angle of the shading - here 30° and 45° - and the direction - the last 
parameter. 
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Line Searching with DERIVE 

George Collie, Dundee, Scotland, g.collie@abertay-dundee.ac.uk 

Line-searching is a basic technique in numerical optimisation. A search direction is obtained say 
\d v d 2 ] when minimising f(v,y) or [ d v d v d 3 ] when minimising a function of three variables. 


Searching along this direction from a starting point at (jc 0 , y 0 ) { or (jc 0 , y 0 , z 0 ) } is carried out by 
minimising f(x 0 + t-d v y 0 + t-d 2 ), {or f(x 0 + t-d v y 0 + t-d v z 0 + t-d 3 )} which is a function of one 
variable t and is called LS(t) in OPT2D_util and OPT3D_util. 


The two simplest methods of obtaining a search direction for minimisation are : 

(1) Steepest descent { - Grad F(jc 0 , y 0 ) or - Grad F(jc 0 , y 0 , z 0 ) } 

(2) The Newton method { - FT 1 Grad F(x 0 , y 0 ) or - FT 1 - Grad F(x 0 , y Q9 z 0 ) } 

The three methods of numerical optimisation for a function of one variable used in OPT2D_util.MTH 
and OPT3D util.mth are: 


(a) Golden section 

(b) Fibonacci search 

(c) Quadratic interpolation 

The steepest descent method ( STEEPD ) uses the direction in which the function has the largest pos- 
sible rate of decrease at (jc 0 , yo ? Zo)- The Newton method ( NEWTD ) uses the Hessian matrix of second 
derivatives to fit a quadratic function to the function and points to the exact minimum of this quadratic 
function. The procedure NEWTSTEP |jc n+1 = x n - H~ l • V/(jc w ) } gives the result of the version of the 

method with no line searches and is equivalent to setting a value of t = 1 in the linesearch LS(t) when 
NEWT2D has been used. The first two univariate search methods are based on finding LS(t) at 2 
points K and L say inside some given range ( A, B ) where A < K < L < B 

Assuming the function is unimodular in [ A, B ]: 
If/(L) </( K) the interval [A, B] becomes 
[K, B]; otherwise it becomes [A, L]. 

When the method is applied to the new range, it 
would be useful if the best value (which is L in 
[K,B] or K in [A,L] could be re-used. 

If AL/AB and KB/AB are each equal to the golden section ratio ( V5 - 1 )/2 ~ 0.618034, then the 
interval is always reduced in the same ratio and the best value can be re-used. 

Similarly, if AL/AB and KB/AB each equal the ratio F n .\JF 2n of two consecutive Fibonacci numbers, 
then this best value can be re-used in the next interval with the next ratio becoming F n _ 2 /F n -i. 

In OPT2D_util.MTH, although the values will be re-used, they will also be re-calculated which would 
not be necessary if the method was programmed in FORTRAN or PASCAL. 

The third method of univariate searching fits a quadratic P t 2 + Q t + R through three points 
\tu LS(/ l i)], \t 2 , LS(£ 2 )], \h, LS(/ l 3 )] and calculates the value t* = - Q/(2P) where the turning point of the 
quadratic lies. The method cannot be guaranteed to produce a better value than any of the t\ but 
(especially if t\<t 2 < t 2 and LS(£ 2 ) < LS(f,), LS(/ 2 ) < LS(^) ), the method often produces a good value. 
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The session printout will show that GOLDEN and FIBON include some hints (programmed using 
IF-statements) for choosing a new range of parameter values, although in practice the best technique is 
usually to accept the first reasonable improvement found). 

OPT3D_util.MTH is an obvious extension of OPT2D_util.MTH with most of the procedures having 
to be changed, although FIBON, GOLDEN and FIBONACCI are unchanged. 


#1 

#2 

#3 

#4 

#5 

#6 


A set of procedures to enable line-searching to minimise F(x,y) 
[F(x, y) x_0 y_0 :=] 

SETX0(x0, yO) := [x_0 := xO, y_0 := yO , Function Value = , F(x0, yO)] 
SETXO(x0, yO) := [x_0 := xO, y_Q := yO, Function Value = , F(x0, yO)] 




tl tZ tE 

#12: QUADSEARCHftl, t2, tB) := t tl tZ tE t4 = t_s 

. LS(t) LS(tl) LS(t2) LS(tB) 

New t 

2 

LSCt2).(tl - tl ■ (t2 + tE) + t2.tB) + LS(tE) ■ (tE - tl) ■ (tl - t2) 

2 ■ (LS(tl) ■ (t2 - tE) + LS(t2).(tE - tl) + LS(tE).(tl - t2)) 

LS(t4) = LS(t_s) 

S - 1 

#13: ga := 
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a K 

t a kk := b - ga-(b - a) 

#14: GOLDENCa, b) := 

LS(t) LS(a) LS(kk) 

IFCLSCkk) < LSOD, New a = a, ) IFCLSCkk) < LS(ll), Best point?, New a = 

L b 

11 :z a + ga- (b - a) b 

LSC11) LS(b) 

kk) IFCLSCkk) < LS(ll), New b =11, Best point?) IF(LS(kk) < LS(ll), , New b = b) 

#15: fb := [1, 1, 2, 3, 5, 3, 13, 21, 34, 55, 39, 144, 233, 377, 610, 937, 1597, 2534] 

FIB O') := fb 
#16: j 



K 

FIBCn - 2>(b - a) 

kk :: a + 

FIBCn) 

LS(kk) 

IF(LS(kk) < LS(ll), Best point?, New a = 


L b 

FIBCn - 1)-Cb - a) 

11 := a + b 

FIBCn) 

LSC11) LSCb) 

kk) IFCLSCkk) < LSC11), New b =11, Best point?) IFCLSCkk) < LSC11), , New b = b) 


FCx,y) = FCx, y) 

#13: check := (xO, y 0) = [x_0, y_0] 

Direction = [d_l, d_2] 


The example is following now: 

#1 : LOADC I : \D0KUS\DNLs\DNL97\MTH26\0PT2D_ut i 1 . mth) 

4 4 

#2: F(x, y) := x + 2-x + y 

#3: SETX0C1, 1) = [x_G := 1, y_0 = 1, Function Value =, 4] 

#4: steepd = [d_l := -6, d_2 = -4] 


#5: LS(1) = 696 

#6: QUADSEARCHCO, 0.5, 1) 




tl 

t2 

t3 


New t 

#7: 

t 

0 

0.5 

1 

t4 : 

= t_s := 0.2433234421 


. LS( t) 

4 

i — ■ 

UJ 

696 

LS(t4 

) = -0.3751293339 

#3: ( 

GETNEWPTt 

) 

;o. 

( 

24332 

34421) 

Y 

FCX, Y) 

#9: 

. -0.459^ 

HO 

6528 

0.0267 

0623145 

-0.3751293339 . 

#10: ' 

5ETX0( -0 . 

,45 

99406 

523, 0. 

0267062 

3145) 
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#11: 

[x_0 -0 

.4599406528, y_0 

:= 0.02670623145, Function Value 

#12: 

steepd 



#13: 

[d_l := -1 

.610806674, d_2 : 

" 5 1 

:= - 7.618997253-10 J 

#14: 

LSC0.25) 

= -1.171522271 


#15: 

L3C0.5) = 

0.03281904911 


#16: 

GOLDENC0, 

0.5) 


#17: 

GETNEWPT ( 

0.1909830056) = 

X Y 

_ -0.767577353 0.02669168046 -1 

#18: 

5ETX0C-0. 

767577353, 0.02669163046) 

#19: 

[x_0 := -0 

.767577353, y_0 

: = 0.02 6691 68046, Fun c t i on Va 1 u e 

#20: 

steepd = 

[dJL := -0.1910504741, d_2 := -0.00007506550317] 

#21: 

LSC0.2) = 

-1.189992455 


#22: 

FIB0NC5 , 

o 

o 

UJ 

w 



#23: 


F(X,Y) 


t 

LSCO 


a K L b 

0 kk := 0.12 11 := 0.18 0.3 

-1 . 188027026 -1 . 19051175 -1 . 190290231 -1 . 186775 509 

New a = 0 Best point 7 New b = 0.18 


In DERIVE 6 we can plot the function in the 
3D-plot window and investigate the minimum 
applying the trace mode. 

The intersection point of the parameter 
curves is given with (-0.8, 0, -1.904), which 
is pretty close to the minimum from above. 



Start and end of session2.dfw which searches for the mimimum of F{x,y,z ): 


#1 : LOAD Cl : \DOKUS\DNLs\ DNL97\MTH26\0PT3 D_ut i 1 . mth) 

#2: Precision := Approximate 

2 2 4 

#3: F(x, y, z) := x + 2-y + z + 3-x-y + 2-z 

#4: 5ETX0(1, 1, 2) 

#5: [x_0 := 1, y_0 := 1, z_0 := 2, Function Value = , 26] 

r -i3 -i3 

#19: |_x_0 := - 2.85874457*10 , y_0 := - 2 . 85874457 -10 , z_0 := 0 . 6176279859 , 

Function Value = , 1. 38077100$] 
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" Special Shaded Planes" 

Sebastiano Cappuccio, Forli, Italy, scappucc@spfo.unibo.it 
and Nurit Zehavi, Israel, ntzehavi@wiccmail.weizmann.ac.il 


On 3 April I received an e-mail from Sebastiano Cappuccio, Italy which extended to a fruitful 
discussion within the DERIVE Newsgroup involving other DERIVIANS, especially Nurit 
Zehavi from Israel. Let's begin the story with Sebastiano's email which might be of interest 
for many of you. I must admit, that the "trick with the grid" was absolutely new for me - a 
shame, indeed. Josef 

With DERIVE n(n> 3), if you plot an univariate or bivariate identity (e.g. x = x, or x + y = x + y), you 
get a ’’shaded plane” with a ’’grid” depending on the Accuracy [1] . 

It is easy to explain this ( DERIVE scans the whole plane in plotting implicit equations). If you want to 
plot the ’’grid” only for positive x, you have to plot an identity true only for positive x: 

abs (x) = x. 

Thus, if you want a ’’shaded” area between the straight lines x = 2 and x = -2, you can plot: 

abs ( (2-x) (x+2)) = (2-x) (x+2). 

If you want to get the ’’grid” between two curves u and v, you have to plot: 

abs ( (u-y) (y-v) ) = (u-y) (y-v) . 

For example try with 

U:=sin(x), v:= 0, orwith u:=x A 2-4, v:=0. 

My problem is: how can I get a shaded area between two curves AND within an interval (for example 
between x 2 - 4 and 0, but only with -2 < x < 2? It seems that DERIVE does not accept logical opera- 
tors in plotting identities. Can someone help me? Thank you. 

As you might imagine, I was very interested to see how this works and probably to find a 
solution for Sebastiano 's problem. See some screen shots: 



Screen shots from DERIVE for DOS 


I don’t know the reason for some "big” meshes appearing in the grid. 

[1] You cannot set the Accuracy in DERIVE 6, so you cannot influence the densitiy of the 
grid! One DERIVE 6 plot will follow after the next DfD plot from 1997, Josef 
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More screen shots from DERIVE for DOS 



Examples 4 reproduced with DERIVE 6 


Some days later we (Sebastiano and I) found an answer from Nurit in our electronic mailbox: 


Hello Sebastiano, 

There is a solution to your problem using DERIVE: 


a : = x rt 2 - pi A 2 
set! : =0=abs (a) +a 
u : =SIN (x) 
v: =0 

b : = (u-y) * (y-v) 
set2 : =0=abs (b) -b 

plot setl 
plot set2 

now simplify 

setl + set2 

and plot the intersection 
of setl and set2. 



All the best, 

Nurit Zehavi and Giora Mann 
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By the way I have just sent an abstract to David Sjostrand in Sweden for the “FUN”-Conference, 
entitled 


FUNtastic Implicit Plotting. 

We (Giora Mann is the co-author have been using games on dots, lines and half-planes aappyling 
expressions similar to the examples presented by Sebastiano Cappuccio. 


Nurit’s next mail arrived on 18 May 1997: 

Several DERIVE - news members mentioned that they couldn't replicate my example and asked me to 
elaborate. 

I will illustrate the method by using another 'nice' example (taken form our teacher textbook): 


#1: Define two circles 

2 

#2: circl := (Ex, y] - [-1, 2]) - 9 

2 

#3: circ2 := (Ex, y] - El, 2]) - 9 

#4: Plot the next two expressions: 

#5: 0 = | circl | + circl 

#6: 0 = | ci rc2 | + circ2 

#7: To plot the intersection: 

#S: 0 = ( | circl | + circl) + (|circ2| + circ2) 

#9: To plot tiie union: 

#10 : 0 = ( | c i r cl | + c i r cl) - ( | c l r c2 | + c i r c2) 

#11: Predict the results of: 

#12: 0 = C | circl | - circl) + (|circ2| - circ2) 

#13: 0 = C | circl | - circl) ■( |circ2| - circ2) 



You can see the plots produced with DERIVE 6. The 
original plots from 1997 are on the next page. 
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Nurit wrote: 

Now I have a question to SWHH: If I ask to plot circl (and not 0 = circl), then the line y = 0 is plotted 
without considering for circl; we have to simplify first. WHY? 

23 May 1997: Another mail from Sebastiano flatters into my mailbox. Sebastiano offers 
additional information. 

The idea of the trick for shaded areas started when a student of mine tried to plot the y = x function, 
but for a typing error, he plotted x = x. I think that it occurs because DERIVE scans the whole plane in 
plotting implicit equations and connects points satisfying the condition (that is all the points of the 
plane, according to the Accuracy set): setting a lower Accuracy, you get a ’’large meshed net” Of 
course this tricks works only with DERIVE 3.x, because former releases cannot plot implicit graphs. 

You certainly know Nurit's messages, maybe you do not know the following messages with great 
ideas from England and Chile. 

From John Alexiou 

DERIVE has a step function STEP (x) = 1 if x > 0 and 0 otherwise. You can use this as a condition. 
Try STEP (4 - x A 2 - y A 2) = 1 and it will plot the region inside a circle. 

Now if you plot STEP ( 4 - x A 2 - y A 2) * STEP(1 - x A 2) = 1 then it will only shade the 
region with -1 < x < 1 which belongs to the circle with radius r = 2. 

I tried and here are the results (See first the respective DERIVE 6 plots): 
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STEP (4-x A 2-y A 2) =1 
STEP (l-x A 2)=l 
STEP (4-x A 2-y A 2) =0 
STEP (x A 2+y A 2-4) =1 
"The intersection:" 

STEP (4-x A 2-y A 2) *STEP (l-x A 2) =1 
"or" 

STEP (4-x A 2-y A 2) +STEP (l-x A 2) =2 
"The XOR- Operation: " 

STEP (4-x A 2-y A 2) +STEP (l-x A 2) =1 
"Plot the next two expressions in one Plot Window:" 
STEP (4-x A 2-y A 2) -STEP (l-x A 2) =0 
STEP (4-x A 2-y A 2) +STEP (l-x A 2) =1 
"How to produce the union?" 

John’s STEP procedure demonstration 



From Claudio del-Pino O. 

You could try this: 

ABS (SQRT (x* (1 - x) * (u - y) * (y - v) ) = SQRT (x* (1 - x) * (u - y) * (y - v) 

in order to obtain a grid between the curves, the x-axis form 0 to 1 . 1 hope that this will work always. 



Compare DfD (left) and DERIVE 6 (right)!! 


DERIVE 6 provides ready made functions for shading regions and it is easy to produce the 
grids by plotting appropriate inequalities. You can look at this on the next page. 

This was Sebastiano's last message and I was eager to check these suggestions, too. 

As you can see in another article in this DNL I was busy to visualize systems of linear 
inequalities and the tools to shade areas presented here are quite different from mine. 
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Which region in the x-y-piane is described by the following system of linear inequalities? 


y > -2 
x - y > -2 
x -3y < 6 


n I ; [vl ;= -Z - |L rZ := y - * - Z, r-3 ; = k * 3 lj - 61 

artl ft = |rl| - rl 

13: iset2 := ft = |t2| * rZ 

M-: := ft = |r3| * i-3 

86: o - |k - y * 2| * |k * 3 y - fs[| + y * Z| + 3 y - 10 


4 

-4 

i& 12 





Find the region of feasible solutions considering the following restrictions (from Linear 
Programming): 


X + 4 v < 24 
2x - 3 y < 23 
3„\ + 2 y < 27 

y > 0 



Shading areas with DERIVE 6 using the 
implemented function and applying Boolean 
expressions (which is sometimes of more 
didactical value!) 



AreaBetweenCurves(u, v, x, -4, 4, y) 


- 4, 2-SINCx: 


2 2 

2-(x - y - 4)-SIN(x) - y-(x -y-4)<QA-4<x<4 


-4 < x < 4 a (u<y<vvv<y<u) 
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The SIMPLEX METHOD on the 77-92 

Bruce Chaffee, CA : USA 

This article is intended to introduce an implementation of the simplex method on the TI-92. For a 
more complete discussion of the method the reader is referred to INTRODUCTION TO 
MATHEMATICAL PROGRAMMING , by Winston. 

The ''Simplex Method " 1 as presented here is often referred to as the "Big M method 1 '. The problem is 
to maximize an objective function subject to a set of constraints (see example below). The numbers 
on the right hand side of the constraints must be non-negative. 

Without much difficulty : the non-negativity of the right hand side of the constraints can be overcome. 
The method can also be used to solve minimization problems (See example FOUR). 

EXAMPLE ONE: < CONSTRAINTS ONLY 


Maximize P = 

GOX! 

+ 30X2 + 2 QX 3 

ENTER DATA INTO MATRIX AS FOLLOWS 

Subj. 

Sx 1 

+ — 1 X 3 _<_ 48 

8 

6 

1 48 

To: 

4x-| 

+ 2 X 2 + 3 / 2 x 3 < 20 

4 

2 

3/2 20 


2 x 1 

+ 3/2X2 1 /^X 3 ^ S 

2 

3/2 

1/2 8 



x-1.x2.x3-> 0 

-60 

-30 

-20 0 




0 

0 

0 0 


To enter the problem, open a new matrix named MAT rix EXample 1 _ (MATEX 1 ) using "ARPS 6 : NEW, 

DATA * MATRIX, VARIABLE MATEX 1 , 5 ROWS, 4 COLS. Once you get into the matrix editor you can 
(and will) change the number of rows and columns as you please. 

Enter the coefficients of the constraints in exactly the same order as they appear on the page. 

The objective function , on the other hand : is rewritten as -60xi - 30x2 " 20*3 + P = 0 and entered on 
the next-to-last line. The "P :1 term is in its own column, and never changes during the solution proc- 
ess, so it may be omitted. 
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[H 
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0 

EH 

eh 

EH 
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EH 

EH 



UM 




EJtol 


pga 

» 

s 

s 

OH 




n 

[!■ 

[afl 


[]■ 

UM 

UM 






^H 

^H 

^H 

^H 

^H 













H 





The row of zeros in the last line is required in all problems and will be used during the solution of 
problems with or :; = :1 constraints. For this first problem : which has only "<" constraints : the 
zeros will remain unchanged during the solution process. 


To make room for the additional variables required for the solution, we use the rules: 

For each £ constraint, add ONE column (a slack variable) 

For each J> constraint, add TWO columns (a surplus variable and an artificial variable) 

For each = constraint, add ONE column (an artificial variable) 

Because this problem has 3 < constraints, we must add a total of 3x1 = 3 columns. Position the 

cursor in the last column and then use the keystrokes “F 6 , 1 , 3 ]: a total of 3 times to add the 3 col- 

umns. LEAVE THE COLUMNS BLANK. DO NOT ATTEMPT TO ADD IN THE ADDITIONAL 
VARIABLES. THE PROGRAM WILL DO IT AUTOMATICALLY. 
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[FT 


i) ~ 


f3w 

Calc 


Other 


FE T Ffi 
F’rgnlOlClg-ar a- 


simplex(matexl ,3,0, □> 




fiflF ilUTB 


FUNt frrl 


1 | 1 1 a 1 1 1 a 

r 

L;r.= : 

' : E ‘ J " 

r 3 = c 

r ’ 




IpRELiniNflRV TRELERU 



‘8 

6 1 

i 

0 0 

481 


4 

2 3/2 0 

1 0 

2G 


2 

3/2 1/2 0 

0 1 

8 


-60 

'38 -20 0 

0 0 

G 


.R 

□ 0 

□ 

0 0 

G . 



Irmih m mutb run. 155 HUB 


Now use “diamond" : “home" to return to the home screen, where you can enter: sim- 

plex(matex1,3 : 0 : G) and run the program. The parameters ; ‘3,D,G r tell the program that there are: 


3 constraints of the a <_ yt variety, 

0 constraints of the "_> 1: variety, and 
0 constraints of the ,: = :| variety,. 


Press £: ENTER' 1 now to show the successive matrices (tableaux) as the problem is solved. When the 
bounded solution is found : all the variables (basic and non-basic) are displayed. 
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The final tableau will be stored in the variable named £: mout :i (Matrix OUTput). This allows you to 
save this matrix in a matrix variable for future reference : or simply display it on the home screen for 
further examination (See example FOUR). 

EXAMPLE TWO: MIXED CONSTRAINTS 


MAX P = 2x-| + 3X2 

Subj 1x 1 + 1x 2 <. 10 
To -1x 1 + 1x 2 :> 5 

OXhj + 1x 2 = 7 
xi x 2 >0 


NOTICE THE CONSTRAINTS MUST BE ENTERED IN 
THE ORDER >_, =. 

ahis is 'how the program knows which inequality is 
which type) 


Load this into the matrix :L matex2 :i as follows Insert 4 additional columns as before 
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Because there is one j< constraint one constraint and one = constraint we must add 
1x1col + 1x2cols + 1x1col = 4 columns of zero’s in the same manner as before ( Fg, 1,3). 

Return to the home screen and enter: simplex(matex2,1 ,1,1) 












D— N— L#26 


The TI-92 Corner: 


Simplex Method 


P 51 


After you press “ENTER 11 , you will notice there are now ones in the last row of the PRELIMINARY 
TABLEAU. These represent an “M" in the row above. For example, a 1 in the last row, with a 0 just 
above it represents :£ G + M :i in the objective function. 

Press ENTER again to display the INITIAL TABLEAU. In the leftmost column, a “-2 r in the next-to-last 
row, with a “1 W below it represents a £l -2 + M 11 as the coefficient of x-j in the objective function. 
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Continue pressing ENTER to step through to the solution. 

When a bounded solution is found, the program checks to see if there is a !£ 0 3: at the bottom of a col- 
umn of any non-basic variable and alerts you. F there is a positive value in this column, then you 
may pivot and see the alternate optimal solution. Because it is possible for there to be more than one 
such column : you are offered the option of selecting the column you want to explore - but you must 
be sure there is at least one positive value in the column. 

EXAMPLE THREE : MORE THAN ONE SOLUTION 
MAX P = 2x 1 + 2 x 2 
Subj 1x<| + 1x 2 <. 5 
To + 1X2 <. 2 

1x^ + 3x 2 > 9 
x 1 X2J>0 
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Enter the problem as before using a matrix name of your choice, such as :£ matex3 n . 
Remember to enter the constraints in the order , >_, = . 

You may then return to the home screen and run the program with £! simplex(matex3 : 2,1,G) n 
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The first solution found is P = 10 @ {1 ,4). If you enter £l 4 n to the prompt in the left screen above, the 
alternate solution is found P = 10 @ (3,2). 
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EXAMPLE FOUR: A MINIMIZATION PROBLEM 

ORIGINAL PROBLEM 

MINIMIZE C= 48y 1 + 20y 2 + 8y 3 
subj. 8yn| + Ay 2 +■ 

to 6yi + 2y2 + 3/2y 3 30 

+ 3/2y2 + 1/2y3 j> 20 
Vi - V2- Y3 21 0 


MODIFIED PROBLEM 

MAXIMIZE P = -48x^ - 20x2 " ® x 3 

subj. Bx-| + 4x2 + ^ x 3 — ^0 

to 6 xhj + 2 x 2 + 3/2xg 1 30 

+3/2x2 + 1/2x 3 < 20 
x 1 r x 2 , x 3 21 0 


This minimum problem can be solved by finding the maximum of the MODIFIED problem, where the 
objective function of the modified problem is the negative of that of the original problem. The coeffi- 
cients of the constraints remain the same. The MODIFIED problem is solved as follows: 
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The three _> constraints require 2x3 = 6 additional columns. When solved, this yields the solution to 
the MODIFIED problem. P = -280 @(0,10,10). The solution to the ORIGINAL problem is then C = -P 
= +280 @(0,10,10). 
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When this final tableau is viewed on the home screen, the matrix is slightly too big to fit on the screen : 
but if you highlight the matrix, you can scroll left and right to see the entire matrix (tableau). 
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Example FOUR is, by the way, the DUAL of example ONE : which is why the value of the objective 
function (+280) is the same in both cases. Displaying the final matrices (using MOUT} for both prob- 
lems on the home screen leads to a natural discussion of the DUAL THEOREM and SHADOW 
PRICES. 
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For example, in the left screen above : where column one is visible, notice the values in the last two 
rows. 

In column 1 the bottom two values, “24” above a "0’' ; stand for 24 - OxM = 24, 
in column 4 the bottom two values stand for 2+ OxM = 2, and 
in column 6 the bottom two values stand for 6+ OxM = a. 

Compare these values to the solution to EXAMPLE ONE (P = 280 @(2,0,8), see the screen below). 
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Similar remarks can, of course, be made about val- 
ues in the last two rows of the final tableau of exam- 
ple ONE, and the solution of example FOUR. 

Questions related to sensitivity analysis can now be 
explored. 

How much does the objective function change if ... ? 
How much can a value change and still have the 
same basis remain optimal? 


The program does not replace sensitivity analysis, but hopefully it will make it more understandable. 


It is obvious that there is a tot of programming code copying necessary to have such a powerful toot. I 
do not expect that you are eager to read and to edit 6 full pages of code. So f offer to download the 
program from my school's homepage. And not only that. Linear Programming was and still is on our 
curriculum. The problem is that we often stop after having written down the goal function together with 
the restrictions but we refuse to start that boring transformation algorithm with more variables. Now we 
can use the T! for not too large, but for large enough problems. 

So very soon the idea came up to add a comfortable input routine to avoid working with the Matrix 
Editor, although editing the matnx sometimes will be usefui to make the process clearer. 

Follow my add-on surface for solving the problem below (A mixed Minimumproblem - 
- Blending problem) 
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We call the program simpQ and then we are offered a new menu: 



We choose Option 2: under * * Option 1. under - * and we enter a 1 :New problem under * ■ 

Then we fill in the data beginning with the coefficients of the restrictions and finishing with the coeffi- 
cients of the goal function. Later on we will be asked to save or not to save the problem "mixed". 
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We enter the coefficients of the rows: 
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The initial tableau:, 


and the solution: 
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given part by part: 
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Wo tv we can save the problem and bad and demonstrate/calculate it next time using Option 2: Open in 
menu point ■ Hnput. We only have to know the problem's name (mixed) and fix Maximum or Mini- 
mum . The parameters have j been stored automatically in a list (mixed!). 

If you want to edit the matrix by hand then do it and call Bruce "s SIMPLEX-routine with one more pa- 
rameter which fixes the kind of output For the 4th parameter choose among n, m f d or b for the 4 
options of menu point ■ ? eg SIM PLEX(mixed ,4,3,1 ,n) will immediately return the solution. Josef 


Some open questions for the TI-92 

Some time ago I sent some questions to David Stoutemyer: 

abs(x-5 ) x < 5 simplifies to -(x - 5). but abs(2x-5) x < 5/2 does not simplify? 
it would be useful to switch from within a program to a data-tabie. is that possible? 

it is very useful to lock a program. But it makes no sense to lock a Text file (Script file}. It cannot be 
opened and executed. It would be helpful to lock it against changes , but let it be executed. 

It would be helpful to have one key (or key combination) to have a pause during plotting. Pressing the 
key should continue the plot-ptocedure. In programs I can help myself with a getkeyQ routine , but I 
think a key press would be useful. 


David promised to look into ail the problems. 












